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OPTICAL  SYNTHETIC  APERTURE  TECHNIQUES 


INTRODUCTION 

The  goal  of  this  program  Is  to  develop  optical  synthetic  aperture 
(OSA)  Imaging  concepts  and  techniques. 

Active  synthetic  aperture  radar  systems  have  been  used  at  microwave 
wavelengths  for  many  years  to  produce  quality  Imagery.  Directly 
emulating  such  systems  at  optical  wavelengths  Is  not  adequate.  As  the 
wavelength  of  the  transmitted/received  radiation  decreases,  a  wavelength 
Is  reached  typically  In  the  near  millimeter  wavelength  region,  where  the 
received  signal  can  no  longer  be  directly  detected  and  classically 
demodulated  to  obtain  the  required  phase  Information.  Rather,  It 
becomes  necessary  to  use  heterodyne  and/or  Interferometric  techniques  In 
order  to  obtain  the  phase  Information.  Further,  at  least  presently, 
optical  sources,  such  as  lasers,  cannot  be  waveform  controlled  in  the 
same  manner  as  microwave  sources.  Thus,  different  concepts,  techniques, 
and  hardware  Implementations  are  required  for  OSA  Imaging. 

Ourlng  this  program,  a  number  of  different  aspects  of  OSA  Imaging 
systems  have  been  Investigated  as  described  In  the  next  section  and 
appendices.  They  Include  a  generalization  of  the  theory  describing  SAL 
Imaging,  a  new  transmitter/receiver  configuration,  a  new  optical 
spectrum  analyzer  for  OSA  systems  and  other  applications,  and  the 
development  of  a  specific  hybrid  processing  configuration  for  processing 
OSA  data. 


SUMMARY  OF  RESULTS 


The  theory  for  generating  OSA  Imagery  was  extended  In  Its 
generality  to  Include  arbitrary  postions  of  the  source(s)  and 
recelver(s)  for  active  systems  and  Includes  arbitrary  receiver  positions 
for  passive  systems.  In  general,  given  arbitrary  object  size  and 
motion,  the  required  processing  Is  space  variant.  However,  with  certain 
geometry  and  object  size  limitations,  the  processing  becomes  space 
Invariant.  Far  field  conditions  lead  to  being  able  to  efficiently 
process  the  data  via  Fourier  transform  algorithms.  One  such  case  was 
experimentally  demonstrated  at  a  wavelength  of  1.22  mm.  This  was  the 
first  SAL  Image  produced  at  near  millimeter  wavelengths.  See 
Appendix  A. 

Some  analysis  was  also  done  regarding  "generalized  Imaging"  of 
which  OSA  Imaging  Is  a  subset.  See  Appendix  B. 

Wide  bandwidth  transmitted  signals  are  required  for  high  resolution 
range-doppler  SAL  Imaging.  Mode-locked  lasers  easily  produce  wide 
bandwidth  output  but  the  efficient  detection  and  processing  of  such 
signals  can  be  troublesome.  A  technique  was  developed  to  reduce  the 
required  detection  (receiver)  bandwidth  without  losing  the  high 
resolution  Imaging  capability.  It  consisted  of  not  only  a  mode-locked 
transmitter  laser  but  also  of  a  mode-locked  local  oscillator  laser  In  a 
heterodyne  detection  setup  that  allowed  the  received  signal  to  be 
compressed  Into  a  narrow  bandwidth.  A  bandwidth  reduction  of  over  500 
was  experimentally  demonstrated.  See  Appendix  C. 

A  useful  Instrument  for  obtaining  the  modulation  Information  from  a 
wide  bandwidth  optical  OSA  signal  would  be  a  high  finesse  optical 
spectrum  analyzer  that  spatially  displays  the  temporal  spectrum  with  a 
one-to-one  linear  correspondence  between  the  temporal  frequency  and  the 


spatial  position.  Such  an  output  would  allow  direct  optical  processing 
or  detection  with  a  linear  detector  array  for  digital  processing  without 
requiring  correction  for  a  non-linear  correspondence.  Towards  the 
objective,  a  spectrum  analyzer  was  conceived  and  demonstrated  that  was 
composed  of  a  mcTtl-tapped  optical  fiber  delay  line  In  conjunction  with 
a  holographic  matched  filter.  This  combination  produces  the  desired 
linear  spatial  distribution  of  the  temporal  spectrum  of  the  light  fed 
Into  the  optical  fiber.  See  Appendix  0. 

The  processing  of  the  three  dimensional  signal  history  data  Is,  of 
course,  an  essential  aspect  of  OSA  Imaging.  In  general,  the  processing 
can  be  very  computationally  burdensome.  However,  by  properly  picking 
the  geometry  and  signal  modulation,  the  data  takes  bn  a  form  that  Is 
conveniently  processed  using  backprojectlon  techniques  similar  to  that 
used  with  tomographic  Imaging  systems.  In  particular,  a  hybrid 
(optical-digital)  processor  concept  that  allowed  such  data  to  be 
efficiently  processed  was  initially  conceived  under  this  contract  and 
then  was  experimentally  demonstrated  with  other  support.  The  further 
development  of  such  a  processor  has  the  potential  to  lead  to  a  real-time 
processor  of  OSA  Imagery.  See  Appendix  E.  Also,  the  grating 
Interferometer  type  of  processor  appears  attractive  for  a  similar  type 
of  OSA  processing.  The  advantage  of  the  grating  Interferometer  lies  in 
Its  achromatic  charqcterlcs  which  considerably  reduce  the  requirements 
on  the  source  coherence.  See  Appendix  F. 

CONCLUSION 

The  concepts,  analysis,  and  experimental  results  developed  in  this 
program  have  helped  establish  the  feasibility  of  optical  synthetic 
aperture  Imaging  and  hopefully  has  laid  the  foundation  for  an 
operational  optical  synthetic  aperture  system. 


APPENDIX  A 


SYNTHETIC-APERTURE  TECHNIQUES  AT  NEAR-HILLIMETER  WAVELENGTHS* 

Carl  C.  Aleksoff 

Abstract 


Synthetic  aperture  (SA)  Imaging  techniques  are  considered  for 
Identifying  targets  at  ranges  of  a  few  kilometers  using  near-millimeter 
waves.  A  generalized  concept  of  SA  Imaging  Is  presented.  A  hybrid  SA 
array  technique  Is  proposed  that  would  significantly  diminish  the  motion 
required  for  the  SA  data  generation.  A  scanning  SA  technique  is 
presented  and  analyzed  In  detail.  Images  are  shown  that  were  generated 
using  the  scanning  SA  technique  with  a  near-millimeter  wave  homodyne 
system  operating  at  1.22  mm. 
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Synthetic-aperture  techniques  at  near-millimeter  wavelength* 


Cert  C.  AMueff 

Radar  and  Optics  Division.  Environmental  Research  Institute  oi  Michigan 
p. 0.  Box  8618.  Ann  Arbor.  Michigan  48107 


Abstract 

Synthetic  aperture  (S4)  imaging  techniques  are  considered  for  identifying  targets  at 
ranqes  of  a  few  kilometers  using  near-millimeter  waves.  A  generalized  concept  of  SA  imag¬ 
ine  is  presented.  A  hybrid  SA  array  technique  is  proposed  that  would  significantly  dimin¬ 
ish  the  motion  reouired  for  the  SA  data  generation.  A  scanning  SA  technique  is  presented 
and  analyzed  in  detail.  Images  are  shown  that  were  generated  using  the  scanning  SA  tech¬ 
nique  with  a  near-millimeter  wave, ■ homodyne  system  operating  at  1.22  mm. 

Introduction 

Some  previous  publications  have  indicated  that  near-millimeter  (NMM)  wave  systems  might 
be  applicable  for  imaging  through  fog,  smoke,  or  rain  in  battlefield  situations  where  visi¬ 
ble  or  infrared  wavelength  imaging  systems  are  blinded.!*3  The  goal  is  that  the  NMM  wave 
imaging  system  would  be  of  a  small  enough  size  that  it  can  be  conveniently  mounted  on  a 
vehicle  and  still  have  a  large  enough  aperture  to  enable  recognition  of  tactical  targets 
at  ranges,  up  to  several  kilometers.  According  to  Kruse*,  recognition  of  tactical  vehi¬ 
cles  requires  about  a  0.5  m  resolution  capability. 

The  required  aperture  width  wj  (Rayleigh  criteria)  of  a  classical  one-dimensional 
imagn.g  system  to  achieve  a  resolution  of  p  at  a  range  R  with  wavelength  \  is  given  by 


or  if  the  system  is  a  focused  scanning  type  that  uses  a  common  aperture  to  both  transmit 
and  receive  the  signal,  then  the'aperture  width  is  given  by  haif  of  Eq.  (1)  or 

*2  "2?  < 


Tyoicallv,  the  nmm  wavelenqr 
aoerture  widths  reauired  '•'or 
mm,  and  3  mm  wavelengtns  as 
figure  1  that  operation  cf  a 
ranges  would  reauire  4n  aper 
antenna  (oish)  that  woulc  oe 
tend  to  crive  the  operating 
the  other  end  of'  the  nmm  wav 
reauirement,  which  obviously 
on  a  vehicle,  especially  in 
t°rnative  techniques  suen  as 


hs  are  taken  to  run  from  0.3  mm  to  3  mm.  Figure  1  shows  the 
a  0.5  m  resolution  as  a  function  of  range  f or  0.3  mm,  1.22 
described  by  the  more  optimistic  Ed.  (2).  It  is  clear  from 
scanning  system  at  a  wavc'sngth  of  0.3  mm  and  up  to  3  km 
ture  of  about  1  m,  which  is  about  the  maximum  size  classical 
tolerated.  However,  atmospheric  attenuation  considerations 
wavelength  selection  to  longer  wavelengths3.  Operation  at 
e  region,  i.e.,  3.0  mm  wavelengths,  leads  to  a  10  m  aperture 
implies  an  impractically  sized  classical  antenna,  for  mounting 
battlefield  situations.  Thus,  one  is  lead  to  considering  al- 
syn'chetic  aperture  imaging  systems.  ' 

Synthetic  aperture  techniques 


0y  synthetic  aoerture  (SA)  imaging,  we  imply  that  a  series  of  steps  of  the  type  Illu¬ 
strated  in  Figure  2  are  usee  to  produce  the  imagery,  we  start  with  an  object  (in  ooject 
space)  illuminated  with  a  soecific  wavefront.  The  signal  .scattered  by  the  object  depends 
on  tne  geometrical  oositions  of  the  transmitters,  detectors,  and  object  as  well  as  specific 
scattering  characteristics  of  the  object  ano  the  spectral  prooert ies  of  the  iilumination. 
The  SA  technique  is  deoenoe-t  m  setting  up  the  object  space  such  that  the  scattered  signal 
has  temporal  variations  cue  to  charging  parameters  in  object  space,  wnich  could  be  caused 
by 

1.  Temporal  modulations  i->  the  transmitted  signal  (e.g.,  cnirping  or  pulsing). 

2.  Transmitter  anj/or  detector  motion  (including  switching  between  sources  and/or  de¬ 
tectors  in  arrays). 

3.  Cb fect  motion  (e.g.,  rotating  object). 

u.  Seam  scanning  (e.a. ,  seaming  a  transmitted  wavefront  with  a  rotating  mirror). 

The  scutered  signal  In  conjunction  with  an  arbitrary  modifying  signal  gives  an  electrical 
signal  which  we  will  denote  'a-  ;ie  SA  temporal  signal.  The  modifying  signal  could  be  a 
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local  reference  wave  heterodyned  with  the  mattered  signal  or  some  more  complicated  com¬ 
bination  of  signals,  electrical  and  electromagnetic.  The  one-aimens  ional  temporal  signal 
is ' then  distributed  into'a  multi-dimensional  spatial  space  (scatiai  signal  space)  via'a 
specific  formatting  procecure.  This  spatial  signal  is  what  is  oTten  referred  to  as  the, 
synthetic  aperture  data  and  the  extc.it  of  the  data  defines  the  synthetic  aperture.  The 
spatial  signal  space  might  also  be  referred  to  as  a  hologram,  especially  when  the  signal 
is  i  ecorded  cnto  photographic  emulsion.  However,  the  soatial  signal  coulo  as  well  be  data 
distributed  in  a  cnmputer  memory  where  the  spatial  coordinates  are  jusi  conceptual  aias 
for  proper  data  retrieval  in  the  next  step.  This  next,  and  last  step,  is  processing  the 
spatial  signal  data  to  give  an  image.  The  processing  could  be  optically  implemented 
Fourier  transform  or  a  digitally  implemented  FFT,  holographic  processing,  or  more  compli¬ 
cated  procedures. 

Often  two  different  types  of  object  space  parameter?  can  be  changed  simultaneously  to 
generate  the  data  for  two  dimensional  imaging.  For  example,  the  classical  side-looking  Sfl 
radar  uses  the  combination  of  range  information  derived  by  transmitting  a  temporally  modu¬ 
lated  waveform  and  Doppler  information  derived  from  the  motion  of  the  aircraft.4*6 

Such  a  "side-looking"  SA  system  also  appears  attractive  using  nmm  waves  in  battlefield 
scenarios.  However,  the  actual  type  of  system  is  more  likely  dictated  by  the  types  of  NMM 
wave  sauries,  detectors,  and  components  available  rather  than  any  idealized. consiaerations. 
If  highly  coherent  (phase  stable)  sources  of  adequate  power  output  and  efficiency  (small 
size)  are  available,  then  heterodyne  systems  are  attractive.  If  detectors  and  components 
(e.g.,  mixers)  are  also  available,  then  direct  detection  and  subsequent  mixing  of  the  elec¬ 
trical  signal  could  be  used,  much  as  for  the  microwave  SA  radars.  If  coherence  of  the 
sources  Is  a  problem  or  if  atmospheric  turbulence  is  a  orobrem,  then  interferometric  tech¬ 
niques  could  be  used,  much  as  has  been  demonstrated  at  optical  wevelengths  or  at  radio 
wavelengths  for  radio  astronomy.  The  distinction  between  microwave  (thermal  noise  limited 
systems)  and  optical  (photon  noise  limited  systems)  becomes  blurred  for  nmm  waves  ana  a 
hybrid  system  will  probably  be  required  in  order  to  blend  the  best  of  microwave  and  optical 
techniques. 

Let  us  consider  how  various  SA  techniques  might  be  implemented.  Given  tactical  situa¬ 
tions,  consistent  movement  of  a  target  is  not  likely  and  hence  target  motion  will  not  be 
considered.  However,  modulation  and  source/detector  motion  will  be  considereo. 

In  a  modulated  radar  system,  we  know  that  the  range  .resolution  capability  is  given  by 
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where  c  is  the  free  space  soeed  of  light  and  8  is  the  transmitted  bandwioth.  Thus,  for  a 
0.5  m  resolution  at  a  .1.2  mm  wavelength,  a  bandwidth  of  300  MHz  is  reauire’d,  which  is  only 
a  0.1%  bandwidth  to  carrier  frequency  ratio.  This  bandwidth  appears  to  be  easily  obtaina¬ 
ble  as  a  freauency  chirp  generated  with  a  backward  wave  oscillator  or  as  a  pulse  train  gen¬ 
erated  with  a  mode  locked  laser.  Thus,  ranging  appears  to  be  a  oooo  candidate  for  gather¬ 
ing  SA  data  with  nmm  waves. 


.  As'  with  microwave  SA  radar  systems  operated  from  aircraft,  we  can  also  envision  the 
cross-range  data  bOing  generated  via  the  motion  of  the  vehicle  carrying  the  nmm  wave  SA 
system.  That  is,  Doppler  or  interferometric  information  is  gathered  to  synthesize  the 
aperture  data  In  s  dimension  orthogonal  tu  the  range  direction,  thus  allowing  the  genera¬ 
tion  of  a  two  dimensional  Image.  This  can  be  accomplished  as  illustrated'  in  Figure  3  for 
a  homodyne  (Doppler)  type  of  SA  imaging  system.  A  diverging  NMM  wavefront  Emanates  from  a 
small  aperture  (antenna  or  lens)  and  Illuminates  the  object.  The  scattered  signal  is  re¬ 
ceived  with  the  same  aperture  and  heterodyned.  For  such  a  system  to  straightforwardly  pro¬ 
cess  tne  signal,  the  position  of  the  receiver/transmitter  phase  center  must  typically^e 
known  to  an  accuracy  on  the  order  of  a  tenth  of  the  operating  wavelength  during  its  coher¬ 
ent  integration  time.  This  is  known  as  motion  compenstion  data.  Thus,  for  a  1  mm  SA  sys¬ 
tem,  the  position  knowledge  must  be  on  the  order  of  0.1  mm  for  up  to  10  m  of  travel  This 

type  of  accuracy  can  be  obtained  if  some  reference  scattered  (e.g.,  retro-reflectors)  arc 
*'3dfesf*<3  vehicle  motion  with  an  optical  or  E4M  distance  measuring 

fact’  lt  should  even  be  Possible  to  use  the  nmm  wave  SA  system  itself  to  qene- 
rate  motion  compensation  signals  from  reference  reflectors  that  are  illuminated  at  the 
same  time  as  the  target. 


Figure  A  illustrates  an  alternative  hybrid  SA  array  technique  where  a  long  thinned  array 
(one-dimensional)  Is  used  along  with  physical  motion  to  fill  in  the  gap  between  antennas. 
For  example,  if  9  small  antennas  were  spaced  1  m  apart,  then  the  entire  array  could  be 
moved  slowly  in  a  track  by  1  m.  while  quickly  switching  (electronically)  the  transmission/ 
receiving  across  the  array  to  obtain  a.  filled  10  meters  worth  of  SA  data.''  \ 
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There  are  other  hypric  SA  techniques  that  can  he  usea,  out  as  statec  earlier,  the  most 
effective  tec,'r'iaue  may  end  up  being  more  determined  Dy  hardware  constraints  then  concep¬ 
tual  constraints.7  Experimental  efforts  are  reauiieo  to  determine  the  effectiveness  of 
the  various  SA  concents  considering  the  hardware  constraints,  atmospheric  turbulence  anc 
scattering  effects,  and  ground  clutter  problems.  In  the  next  section,  we  will  Describe  a 
SA  experiment  that  is  in  a  very  early  stage  of  demonstrating  NMM  wave  SA  imaging 
capab il i ties. 


Scanned  SA  imaging 

Consider  first  the  classical  scanned  imaging  system  illustrated  in  figure  5a.  A  focused 
oeam  is  scanned  past  the  target  using  a  rotating  mirror.  The  amplitude  or  intensity  of 
the  received  signal  scattered  from  the  target  is  displayed  on  a  screen  with  its  position 
proportional  to  the  scan  angle.  In  figure  5b,  the  SA  scanned  system  is  schematically  illu¬ 
strated.  However,  in  this  system,  a  diverging  wavefront  is  scanned  past  the  target  using 
the  rotating  mirror.  The  received  signal  (heterodyne  detection)  is  stored  in  a  memory  un¬ 
til  the  entire  target  is  scanned  ana  then  the  signal  is  processed  to  form  the  image.  Both 
i /stems  have  the  same  resolution  limit,  namely  that  given  by  Eg.  (2)  where  «2  is  the  width 
of  the  scanning  mirror.  The  SA  system  usually  can  have  simpler  optics  since  only  a  diverg¬ 
ing  is  oeam  needed  and  it  does  not  have  to  be  focused  to  a  specific  range  with  the  optics. 
However,  time  must  be  spent  to  process  the  SA  data  before  the  image  is  formed.  The  SA 
scanning  system  has  the  advantage  that  it  can  refocus  by  reprocessing  the  data  and  hence 
only  one  scan  is  necessary,  however,  it  must  keeo  the  SA  data  in  memory. 

Let  us  now  develop  the  processing  concept  for  scanned  SA  imaging  (refer  to  figure  6). 

In  this  figure,  S  is  a  point  source^  0  a  point  detector,  Pn  is  the  nth  point  scatterer 
of  N  scatterers  making  up  the  target,  and  M  is  the  scanning  mirror,  from  the  vantage  posi¬ 
tion  of  the  source  and  detector,  the  object  point  appears  to  move  in  a  circle  of  radius  cn 
centered  at  the  mirror  rotation  center.  In  fact,  it  would  move  on  a  sphere  if  the  mirror 
were  a  two-dimensional  scanning  type,  however,  for  simplicity,  we  will  confine  ourselves 
to  only  considering  a  one-dimensional  scan.  In  regards  to  figure  1,  it  is  the  rotation  of 
the  mirror  that  produces  the  "change"  indicated  after  the  first  box.  With  respect  to 
figure  6,  let  the  transmitted  signal  from  S  be  a  diverging  spherical  wave  of  the  form 


(A) 

2tt/i  s  u/c.  The  signal  scattered  back  to  the  detector  is  of  the  form 


exp  { i[wt  -  k(rp  ♦  dn)J}  (5) 


where  a n  is  the  scattering  coefficient  for  the  nth  scatterer,  rn  is  the  distance 
from  che  source  to  the  nth  scatterer,  and  dn  is  the  distance  from  the  detector  to  the 
nth  scattere.'  After  homodyning  with  the  local  oscillator  exp  (iwt),  the  detected  signal 
is  proportional  1? 

N  •’ 

U(t)  =  £on  cos  [2s  (rn  -dn)]  (6) 

n 

where  «n  not  only  accounts  for  the  scattering  efficiency,  but  also  losses  in  the  system 
as  well  as  (rndn)-^  which  we  assume  is  so  slowly  cnanging  that  it  can  be  taken  as  a 
constant.  Also,  on  is  assumed  to  be  real  for  simplicity  siake.  Of  course,  rn  and 
dn  are  functions  of  time  since  they  vary  with  mirror  position.  The  signal  u(t)  is  the 
temporal  signal  referred  to  in  the  second  box  of  figure  1,  while  the  local  oscillator  Is 
the  modifying  signal.  The  next  step  is  to  consider  the  formatting. 

Towards  understanding  the  formatting  process,  we  first  note  from  Eq.  (6)  that  the  terms 
in  the  summation  are  cyclic  in  rn  ♦  dn'.  In  fact,  rn  ♦  dn  *  q  ,  q  =  0,  1,  2,  .  .  ..de¬ 
scribe  a  set  of  ellipsoids  of  revolution  in  space  with  S  and  0  as  foci,  we  will  denote 
these  "imaginary"  surfaces  as  virtual  fringes,  for  a  reason  to  soon  become  evident.  Con¬ 
sider  the  holographic  recording  setup  illustrated  in  figure  7.  The  hologram  is  to  be  re¬ 
corded.  on  the  recording  surface  using  a  spherical  wavefront  converging  to  0  and  a  diverging 
wavefront  emanating  from  S.  The  recording  at  the-  point  P  is  given  by 


-  e^p  [  ioit  -  kr  /  ] 


where,  as  usual,  k  = 
N 


U'(t)  = 


n 


on 

rnd_ 
n  n 
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| A  exp  [-Hut  -  kr)]  ♦  8  exp  ♦  Wr)J| 

»  A2  ♦  82  ♦  2AB  cos  [k(r  ♦  d)]  (7 

Notice  that  the  phase  term  depends  on  r  *  d  and  then  compare  to  Eq.  (6).  The  hulogram  re¬ 
cording  made  under  these  conditions  has  the  same  form  as  the  signal  received  with  the  Sfl 
system.  That  is,  each  term  In  Eq.  (6)  can  be  considered  to  be  a  hologram  recorded  on  a 
circle  of  radius  cn  and  formed  by  the  interference  between  a  spherical  wavefront  converg¬ 
ing  to  the  detector  0  and  a  diverging  wavefront  emanating  from  the  source  S.  Thus,  the, 
concept  of  virtual  fringes  in  the  SA  system  are  derived  from  the  real  fringes  of  this  com¬ 
parison  holographic  recording. 

Keeping  this  holographic  concept  in  mind,  the  signal  u(t)  is  recorded  along  a  circle  in 
recording  space  by  formatting  u(t)  into  the  function  u(p,  a)  via  the  relationship 


U(p,  ♦  )  .  f U(t)6(  p  -  pQ,  a  -  2*m) 


where  ( p,  ♦)  are  polar  coordinates  in  recording  space,  p0  is  an  arbitrary  constant,  <5(p,$ 
is  a  two-dimensional  delta  function,  and  is  the  angular  position  of  the  mirror  at  time 
t  and  hence  is  a  function  of  t.  Recall  that  the  object  will  appear  to  rotate  at  twice  the 
angle  of  the  mirror  rotation  angle.  If  the  mirror  rotates  at  a  constant  angular  rate  of 
fl,  then  we  can  write  that  a*  *  Qt. 

The  processing  step  is. a  holographic  type  of  reconstruction  as  illustrated  in  Tigure  a. 
within  the  limitations  of  the  Fraunhoffer  approximation,  the  image  complex  amplitude  is 
described  oya 


Uj(R,0)  « U(  o,  ♦  )Uc(p,  o  )w(p ,  o  ,  R,0)| 


where  represents  a  two-dimensional  Fourier  transform,  uc(p,  a  )  is  the  reconstruction 

wavefront,  w(  p,  a  ,  R,8)  is  the  propagation  kernel,  and  (R,  9)  are  the  polar  coordinates 
in  image  space,  when  Eq.  (?)  is  substituted  Into  Eq.  (9)  we  obtain 


ut(R,  9).  Sn(R,  9)  ♦  Sn*(R ,  *] 


where 


Sn(R’*  )  *  *3  |V*  *  ?fex(>[¥  (rn  -  V]6*  9  '"o'  ♦  *  20n)dt  [ 


we,  note  that  Sn  is  the  spread  function  corresponding  to  the  object  point  Pf,.  The  qual¬ 
ity  of  the  image  can  oe  determined  from  the  characteristics  of  the  sp^aJ  function.  Aber¬ 
rations,,  for  example,  can  be  analyzed  by  comparing  the  position  of  each  spread  function  to 
that  of  its  corresponding  object  scatters?  point.  Rather  than  expand  and  analyze  Eq.  (ii) 
in  detail  for  arbitrary  configurations,  we  will  consider  a  special  case  that  not  only  is 
straightforward  and  demonstrates  the  concept,  but  corresponds  to  our  experimental  setup. 

As  indicated  in  Figure  9,  consider  the  case  where  the  source  and  detector  are  in  the 
same  position  and  the  wavefronts  collimated.  Hence,  the  virtuai  fringes  have  degenerated 
to  equally  spaced  straight  lines  of  porto <i  1/2.  If  the  mirror  is  rotating  at  constant 
angular  frequency  n,  then  the  received  signal,  as  described  by  Eq.  (4),  becomes 

N 

u(t)  .  on  cos  cfl  cos  (2Slt  -  s,,)  *  o]  (12, 

n  ■ 

where  «  is  a  constant  describing  the  phase  delay  within  the  system  optics.  Note  thet  <t» 
and  sn  are  both  constants.  The  formatting  orocedure  gives  the  function 


u(p ,  a) 


°n  co*  "To '  co*  <  ♦ 


*  •  „)  *  «  A(p, 


where  A(p , a  )  describes  the  effective  aperture  (a  section  of  .circular  arc) 
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A  ( P ,  «  ) 
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and  T  is  the  scan  period.  Notice  that  the  nth  term  in  the  summation  in  Eq.  (13)  describes  a 
sinusoidal  grating  with  amplitude  an,  frequency  2cn/Xpo,  anyular  orientation  8 n,  and  offset 
u.  The  processing  now  degenerates  to  using  uc  =  1  and  *  *  1  in  Eq.  (9)  to  give 

Uj(R,  9 )  =  a(R,  9)  * 


N 

2£|y<*  -  *"•  ®-  8n)  e-ia 
0 

*  dn»(R  -  %.  0  ♦  ®n>  eia] 

where  spread  function  a(R,0)  is  the  Fourier  transform  of  fl(p,  ♦  )  and 


(15) 


(16) 


Eq.  (15)  indicates  that  we  get  the  usual  dual  images  from  the  holographic  process.  Since 
the  center  of  rotation  is  outside  the  object  and  the  object  would  typically  subtend  less 
than  180°,  then  the  two  images  do  not  overlap.  This  is  the  desired  imaging  result. 

Experimental  Results 

Figure  10  schematically  illustrates  the  scanning  nmm  wave  SA  system  that , was  used  to 
generate  our  experimental  data.  This  system  is  part  or  a  MICOM  nmm  wave  radar.10  The 
source  is  a  C13  h3F  wavaguide  laser  pumped  by  a  CO2  laser.  The -expanded  and  collimated  C# 
output  at  1.22  mm  is  scanned  with  a  flat  12  inch  diameter  mirror.  Part  of  the  beam  is 
tapped  off  with  a  beamsplitter,  and  forms  the  local  oscillator  for  homodyne  operation.  The 
return  signal  from  the  target  ilong  with  the  local  oscillator  signal  are  detected  (homo- 
dyned)  at  a  C-oloy  cell  which  has  about  a  OC  to  15  Nr  bandwidth. 


Figure  11(a)  shows  a  plot  of  a  received  signal  from  a  60  mm  diameter  corner-reflector 
that  was  117  cm  from  the  rotation  center  of  the  scanning  mirror.  The  total  scan  angle  was 
17.5°  and  the  signal  lata  was  digitized  at  350  samples  corresponding  to  0.05®  intervals. 

The  reflector  was  within  m«  same  room  as  the  nmm  wave  system  in  order  to  avoid  atmospheric 
degradations,  at  lea't  in  the  initial  experiments.  Compare  the  signal  in  Figure  11(a)  to 
that  in  Figure  12(a)  which  is  a  computer  simulated  signal  from  an  idealized  point  scatterer 
illuminated  with  infinite  wictn  plane  waves.  The  actual  return  is  obviously  not  * ideal *. 
me  non-symmet ric  effects  may  ^ue  to  polarization  effects  from  the  comer-reflector 
since  the  reflector  was  not  symmetrically  orientated  in  the  scanning  direction.  Also, 
feedback  effects  from  the  reflector  into  the  laser  were  evident.  Further,  the  output  beam 
was  somewhat  Gaussian  shaped  with  a  full  width  at  half  power  points  of  about  ICO  mm. 


For  our  computations  the  discrete  version  of  Eg. 
namely 


Uj(R ,  9) 


m«  1 


-1 
0  ’ 


mi.C  ) 


cos 


cos(  9  - 


(9)  was  used  for  unity  magnification, 


where  R0  *  1170  mm  and  4$ «  0.05°,  ana  m  <  3 5 i .  The  Image  snpwn.  in  Figure  11(d)  was  ob¬ 
tained  by  taking  a  cut  through  image  space  along  an  arc  of  radius  R0  as  described  by  the 
image  samples 

!q  *  I W  <3A*)I  2  a  *  1 . 151 


Figure  12(b)  is  a  construction  of  me  image  formed  from  the  "Ideal"  signal.  In  comparison, 
the  experimental  results  snow  ime  expected  narrow  spike  from  the  corner  reflector  out  with 
a  nigh  level  of  background  noise.  :t  is  expected  that  some  additional  processing,  such  as 
bandpass  filtering  the  signal,  xnould  improve  the  quality  of  the  image.  However,  it  snould 
be  noted,  that  the  theoretical  r-* solution  of  7  mm  at  tne  target  for  a  100  mm  wloe  Gaussian 
wavefront  was  achieved  as  demonstrated  cy  the  experimental  results. 
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Figure  1.  One  dimensional  aperture  width  required 
for  1/2  meter  resolution  as  a  function  of  range. 


Figure  4.  Hybrid  synthetic  aperture  techniques 
where  a  thin  array  Is  fllled-ln  by  translating 
the  array. 


Simulated  results  for  the  same  parameters  as  for  Figure  11  assuming  an  Ideal  point  scatterer.  (a)  The  temporal  signal. 
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GENERALIZED  IMAGING 
C.C.  Aleksoff 


The  goal  of  this  paptr  Is  to  start  our  thinking  process  about  what 
are  the  fundamental  elements  of  an  Imaging  system  and  what  do  we 
understand  of  these  elements.  The  workshop  goal  Is  to  Identify  those 
areas  where  research  Is  warranted. 

The  word  "Imaging"  usually  brings  to  mind  the  human  type  of  visual 
system  consisting  of  a  lens  and  an  array  of  photoreceptors.  Processing 
Is  then  done  on  the  entire  array  of  photoreceptor  outputs.  However,  It 
Is  clear  that  Insects  and  other  creatures  have  significant  variations  on 
this  theme.  Man  made  scanning  systems  and  synthetic  aperture  systems 
also  do  not  fit  In  that  mold.  It  appears  that  In  order  to  properly 
cover  all  the  possibilities  we  need  to  understand  the  whole  chain  of 
imaging,  from  the  source  of  energy  and  Its  properties  to  the  very 
utilization  of  the  final  data.  What  Is  It  that  has  not  been  exploited? 
It  Is  the  much  broader  class  of  generalized  Imaging  that  we  wish  to 
attack. 

Imaging  can  be  usefully  subdivided  for  analysis  into  the  set  of 
categories  Indicated  by  the  boxes  In  Figure  1.  They  consist  of 

t.  Illumination  (of  the  scene) 

2.  Scene 

3.  Channel  (scene  to  collector  Interval) 


Figure  1 


c.yyryy-S: 


4.  Measurement  system  (which  consists  of  collectors  and 
.  detectors  on  platforms) 

5.  Inversion  (to  generate  an  "Image"  from  the  detected  Signal) 

6.  A  priori  knowledge  (to  help  the  inversion  or  utilization) 

7.  "Image"  utilization  (generation  of  a  useful  output  as  an  Input 
for  a  decision  making  process) 

In  Table  7  these  categories  are  further  subdivided  Into  descriptors  that 
Indicate  Important  attributes  or  functions  of  the  categories.  These 
descriptors  help  guide  our  understanding  of  the  extent  of  each  category. 
(Notice  that  In  Table  I  with  respect  to  Figure  1  that  the  Illumination 
has  been  Incorporated  Into  the  scene  as  a  descriptor  and  that  a  priori 
knowledge  has  been  left  out,  but  Is  Intrinsic  with  many  of  the 
algorithms  associated  with  Inversion). 


One  way  of  utilizing  Figure  1  and  Table  I  Is  to  track  a  particular 
system  through  the  various  categories  apd  descriptors.  For  example  let 
us  consider  the  human  visual  system. 


Illuminator: 

Scene: 

Channel : 

Measurement  System: 


Inversion:. 

A  priori  knowledge: 
Utilization: 


Sun  or  other  source  of  visible  light 
Reflectivity  and  movement 
Atmospheric 

Collector;  imaging  lens  (two) 

Detector;  photoreceptive  array  (two) 
Platform;  eyeball  ♦  head  *  body  movement 

Direct  ♦  binocular  effects 

Relative  scale  of  objects  for  depth  cue 

Mensuration,  change  (motion)  detection,  etc. 


We  certainly  have  not  considered  all  the  possibilities  In  the  above 
example,  but  It  does  start  us  thinking  In  a  systematic  way  about  the 
system. 
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Let  us  consider  another  example,  namely.  Aperture  Synthetic  Radio 
Astronomy. 


Illumination:  None 

Scene:  Self  emissive,  spatial  Incoherence 


Channel:  Free  space  ♦  atmosphere 

Measurement  System:  Collector;  antenna  array 

Detector;  microwave  sensor 
Platform;  earth  (rotation  essential) 
Inversion:  Fourier  tranform  of  correlation  between 

sensor  outputs  to  produce  a  sky  brightness 
map  (Image) 


A  priori  knowledge:  Non-negative  sky  brightness 
Utilization:  Source  distribution,  mensuration 

Reprocessing:  Strong  source  removal  ("clean*  algorithm) 


In  many  respects  this  system  is  very  different  than  the  human  visual 
system  yet  produces  an  "Image."  It  would  be  Informative  to  track 
through  these  categories  and  descriptors  for  other  creature  Imaging 
systems,  but  I  will  leave  the  tracking  up  to  those  who  are  experts  In 
such  systems. 


As  Indicated  at  the  start,  this  paper  was  Intended  to  stimulate 
thought  about  Imaging  and  not  necessarily  present  answers.  However, 
some  observations  will  be  made  here.  It  1$  clear  that  a  priori 
Information  pi *ys  a  key  role  In  extending  the  capabilities  of  many 
systems  that  gather  Incomplete  data  sets.  Modern  computational 
capabilities  often  allow  this  a  priori  data  to  be  used  where  It  was 
Impractical  before.  An  example  of  this  Is  Iterative  processing  which 
constructs  an  Image  from  only  Fourier  modulus  data  and  a  non-negativity 
constraint.  Modem  computational  capabilities  are  often  essential,  as 
In  the  ease  of  the  Aperture  Synthesis  Radio  Astronomy  imaging  example 


given  above.  Typically,  In  this  case  an  Image  Is  formed  only  after 
Intensive  data  processing  (Inversion). 

Are  there  other  key  elements? 
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APPENDIX  C 


LASER  HETERODYNE  RAN6IN6  USING  A  MODE-LOCKED 
LOCAL-OSCILLATOR  FOR  BANDWIDTH  COMPRESSION* 


Carl  C.  Aleksoff  and  L.  M.  Peterson 


Abstract 


It  Is  shown  that  a  mode-locked  laser  can  be  used  as  a  local 
oscillator  (LO)  In  heterodyne  ranging  to  bandwidth  compress  the  received 
signal.  The  amount  of  bandwidth  compression  depends  on  the  difference 
In  the  Intermode  frequencies  between  the  LO  laser  and  the  transmitter 
laser,  which  can  also  be  a  mode-locked  laser.  Since  the  heterodyning 
occurs  on  the  detector,  It  Is  only  necessary  to  use  a  narrow  bandpasi 
receiver  that  passes  the  compressed  bandwidth  to  obtain  ranging 
equivalent  to  that  obtlned  by  a  wide  bandwidth  receiver  that  passes  the 
full  transmitter  bandwidth.  Bandwidth  compression  of  over  500  times  and 
one  foot  range  resolution  was  exp  rlmentally  demonstrated  with  a  1.5  GHz 
He-Ne  laser  mode-locked  transmitter  signal  was  reduced  to  less  than  3 
HHz  at  detection. 


to  Optical  Engineering  for  publication.  Parts  of  this 


paper  were  also  presented  at  the  Los  Alamos  Conference  on  Optics  '83  and 


published  as  "8andw1dth  Compression  Using  a  Mode-Locked  Local  Oscillator 
Laser"  In  the  Proc.  SPIE,  Vol.  380,  1983,  pp.  239-247. 
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Abstract 

It  has  been  demonstrated  in  a  laboratory  experiment  that  a  mode-locked  laser  can  be 
used  as  a  local  oscillator  'L0)  in  a  heterodvne  ranging  system  to  sample  or  bandwidth  com- 
cress  the  received  signal.  The  bandwidth  compression  deoends  on  tha  difference  in  the 
intermode  frequency  between  the  L0  laser  and  the  transmitter  laser.  Bandwidth  compression 
of  over  500  times  was  demonstrated,  i.e.,  a  1.5  GHz  Ke-Ne  laser  transmitter  signal  was  re¬ 
duced  to  3  MHz  upon  detection.  Since  the  mixing  occurs  on  Che  surface  of  the  detector  only 
a  narrow  bandoass  receiver  that  oasses  the  compressed  bandwidth  is  necessary  to  obtain 
ranging  equivalent  to  that  obtained  by  a  wide  bandwidth  receiver  that  passes  Che  full 
transmitter  bandwidth 

Introduction 

In  laser  ranging  systems  and  in  laser  radar  (LAOAR)  the  resolution  of  the  system  lm- 
oroves  as  the  bandwidth  of 'the  transmitted /received  signal  Increases.  Basically  the  best 
range  resolution  &R  for  a  monostatic  LADAR  system  is  given  byl 

AR  -  vm  (II 


where  c  is  the  soeed  of  light  and  B  is  tha  bandwidth  of  che  eransmlcced/recelved  signal. 

For  a  system  that  transmits  pulses  of  width  T,  for  example,  using  mode-locked  lasers,  then 
3  -  T*l  can  be  substituted  into  Eo.  (1).  However,  in  general,  pulses  are  not  required 
and  ocher  waveforms  such  as  linear  chiro  signals  can  be  used  to  minimize  peak  power  levels. 

Tvo<caliv,  che  larger  che  transmitted  bandwidth  the  larger  Che  detector  and  receiver 
bandwidchs  must  be  to  utilize  che  signal  and  obtain  che  best  resolution.  However,  the 
technique  described  in  this  oaoer  alleviates  the  detector  and,  receiver  bandwidth  require¬ 
ment  bv  caking  advantage  of  che  low  bandwidth  occupancy  of  mode-locked  lasers.  That  is, 
lasers  often  emit  only  in  discrete  longitudinal  modes  with  no  emission  between  che  modes. 
For  a  mode- locked  laser  the'  modes  are  phase  locked  together  to  form  che  pulses  characteriz¬ 
ing  such  lasers*.  Further,  since  the  information  ^modulation)  rate  of  return  from  a  tar¬ 

get  is  often  low  then  this  information  can  be  redistributed  to  lie  between  the  modes  in 
such  a  wav  as  to  reduce  Che  overall  bandwidth  requirements  on  che  detector  and  receiver. 

»  ' 

Basics 

The  technique,  which  we  denote  as  differential  mode-locked  laser  OMLL)  ranging,  uses 

t"o  seosrace  mode-locked  lasers.  One  mode-locked  laser  is  che  source  of  che  transmitted 

e-e-gv  and  Che  ocher  is  used  as  Che  lo-al  oscillator  CLOV  in  a  heterodyne  type  of  setups. 
L«e  s^/cv  be  the  transmitted  signal  and  sq(c)  be  the  L0  signal,  and  S^(f)  and  S(j(f)  their 
sdeccra  iFourier  transforms),  respectively .  The  signals  are  related  by  the  Fourier  trans¬ 
fer  relationship 

r  -2rifc 

Sff)  •'  /  sfOe  dt  (2) 


He  mode-locked  spectral  signals  can  be  written  as 

Vn  *  Z*n4ff  *  "r0  *  f0  *4f>< 


5lff)  -  £ 


\i(l  -  nPl  -  f0). 


which  has  che  time  signal  description 


»,ft)  -  V  b_  exof -2 *i fnF.  *  f-Vcl 
1  n  n  10 
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,,here  f o  i*  an  arbitrary  center  oocical  freouenev,  Fo  and  F^  are  Che  incermode  frequence 
differences  for  the  LO  and  transmitter  signals,  resoectively ,  an  and  b_  are  Che  i'comolex) 
mode  amolicudes  of  the  LO  and  transmitter  signals.  respectively',  and  dr  is  che  LO  offset 
freouenev. 

If  s^fti  is  transmitted  and  scattered  off  some  tercet  then  the  received  signal  sjiO  can 
be  written  as 

s2(t>  -  -  cqi  (6> 

where  c„  is  che  round  trio  delay  time  end  e»  is  the  effective  amolitude  scattering  coeffi¬ 
cient  associated  with  che  cth  scatterer  of  the  target.  In  general,  tn  is  a  function  of  time 
if  the  target  is  moving.  Substitution  of  Eq.  (5">  into  Eq.  (61  gives 


s*(tl  - 


h  (tlb.  e*pf-2irt(fft  ♦  nF,  r  d)  ] 


where 


hn(t>  •  51  c  exp[2*i(f,j  ♦  nF^it  ♦  tdl  ( 

q  h  1 

In  che  above  exoressions,  d  is  a  frequency  shift  chat  accounts  for  the  noainei  Doppler 
shift  of  the  target.  The  received  coamlex  soectrua  is  obtained  by  taking  the  Fourier 
transform  fFT'  of  Eq.  (7',  which  gives 

vfi  -  Iwf  -  f0  -  nri  - <*>  < 

n 

where  Hn(f>  is  che  FT  of  hn(t>.  The  quantity  Hft(f>  describes  che  soectral  screed  of  Che 
nth  mode  due  co  che  carsec.  Tvotcallv,  Che  spread  is  limited,  and  we  assume  that  it  is 
limited  to  a  bandwidth  tf,  i.e., 

Hn(fl  -  0  for  | f |  >  *  (l 

The  detected  signal  for  sufficiently  fast  detector  end  receiver  electronics  is  given  bv 

tft'  •  |  S q ( C 1  Sjfti  I  2 


•  *0*0  *•  *2*2  *  *0*2  *  *2*0  fU' 

where  *  reoresents  comolex  conjugation.  The  detected  spectrum  is  Chen  che  FT  of  Eq.  (ID, 
which  is 


*  *00  +  x22  *  *02  *  x20 


where  che  correlation  soeecra  are  defined  via 


xan  *  -  *.cn  *  Vf> 


After  SMhscltucfon  of  Eq.  (V  and  Eq.  the  various  correlation  spectra  are 


-  r  fidf’ 

\  m  .  n 

i  and  Eq.  the  various  correlation  speed 

n0  "  Z  5*n*C  S'"o  *  "r0  * 

<22  -  II"Xyf  -  nry  *  vf  *  "V 

>•20  •  W  Vk<^f  *  <  *  kF0  •  nFl  .  (1 
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These  cor-elstion  soeccra  are  shown  in  Figure  1  for  a  simple  illustrative  case  of  two  mode- 
locked  lasers  with  just  three  equal  amolitude  modes,  i.e., 

s0(f)  -  £  '  fo  *  nFo  ■  fl81 


-  £  iff  -  f0  -  nfp 
n«-l 


and  a  return  givan  b* 


S2fn  **£  rece  *  f0  *  <1  -  nFj^/Wl  (20) 

n*-l 

The  ooint  to  note  is  that  under  certain  conditions  the  received  spectra  is  available  as  a 
band  comoressed  spectra  within  the  detected  spectra  in  the  region  W/2  <  f  <  Fi/2,  i.e., 
as  illustrated  in  Figure  l(g'  within  the  region  labeled  "bandpass'*  (compare  with  Figure 
Ka)).  To  quantify  this  observation  let 

F1  <  Fo 

H  (f>  -  0  for  If!  >  W/2  (21) 

n 

4F  -  F0  -  Ft 
f  -  Ft/4 

The  bandpass  region  fgp  then  ranges  from  W/2  to  T\/Z.  We  can  gee 

F,  -  W 

*  *  '7LY'~  (22) 

modes  into  the  oassband  without  overlap.  The  limit  is  when  we  let  AF  •  W  then 

F,  -  W 

N  -  (23' 

tf  the  above  conditions  are  satisfied  Chen  we  can  bandpass  Che  terms  in  Eq.  (It*)  for  which 
n  »  k  to  obtain  the  reel  part  of 

h<t'  "  ZbJXX^nAF  -  f)  (24) 

n 

vh»-e  AF  ■  Fg  -  F^  amd  fg  *  Af  -  d.  Comparison  of  this  Equation  with  Eq.  (9)  shows  that 
all  the  desired  terms  are  oreaene,  however  there  Is  a  difference  in  chat  a  weighting  an 
exists  and  che  frequency  direction  is  reversed.  The  conjugation  is  not  significant  yince 
the  terms  are  seoayaced  and  only  che  real  part  is  used.  The  new  carrier  frequency  fn  and 
che  new  tncermode  spacing  AF  are  control  Table  via  Af  and  Fg.  If  the  10  laser  has  uniform 
amolitude  modes,  i.e.,  an  -  l,  then  che  filtered  output  is  exactly  of  the  same  form  as 
silt'  exceoc  time  is  reversed,  which  is  not  significant.,  What  is  imoorcanc  is  chat  che 
signal  mode  spacing  has  been  compressed  by  the  ratio 


which  is  the  desired  result. 


Excar itnencal  Results 


The  experimental  setuo  is  shown  in  Figure  2.  Two  Ke-Ne  lasers  were  mode-locked  to  che 
irequeneies  Fg  and  F|,  which  ior  the  results  oresencsd  were  aboue  73.94  MHx  and  73.80  MHz, 
•esoectivelv .  The  mode  locking  was  accomplished  via  intracavity  acoustic  cells  driven  at 
half  the  mode  lock  frequencv2.  The  mode-lock  difference  frequency  wee  aF  •  Fq  *  Fi. 

■  140  kHz.  The  signal  derived  from  che  mode  lock  *F  oeciltacors  vie  a  mixer  Ye  half  this 


frequency  and  Is  used  co  synchronize  Che  osctlloscooe  connected  co  the  output  detector. 

The  laser  operating  at  frequency  F^  was  used  as  the  transmitter  source  fits  light 
illuminated  che  target! ,  and  the  laser  ooerating  at  frequency  Fq  was  used  as  the  LO  laser 
source.  The  LO  laswr  tied  one  of  Its  cavity  mirrors  mounted  on  a  piezoelectric  translator 
in  order  to  control  the  laser  cevlf»  length  which  determine  the  LO  offset  frequency  Af. 

A  half  wavelength  change  in  cavity  length  allows  a  full  74  MHz  of  change  In  Af.  Both 
lasers  were  monitored  by  a  high  speed  detector  connected  to  sampling  oscilloscopes  and  were 
also  monitored  via  an  optical  spectrum  analyzer. 

Light  scattered  from  the  target  was  combined  with  che  LO  light  by  a  beamsplitter  and 
che  two  resulting  combinations  were  detected  wtth  two  detectors.  One  detector,  a  fast  >2 
GHz  cutoff  frequency  detector,  was  connected  co  an  electrical  spectrvas  analyzer  co  observe 
the  entire  electrical  spectrum,  while  che  ocher  detector,  a  alow  <70  MHz  cutoff  frequency 
detector  wee  fed  to  an  oscilloscope  via  a  filter  to  observe  the  desired  filtered  time  sig¬ 
nal.  It  is  noted  that  this  is  a  heterodyne  technique  and  hence,  that  the  wavefronts  of 
che  two  light  fields  have  co  be  aligned  accurately. 

Figures  3  and  4  show  some  of  che  results.  Figure  3  presents  che  results  for  e  target 
composed  of  s  single  point  scatterer,  and  Figure  4  for  two  point  scatterers  separated  by 
about  14  inches.  Figures  3 (a!  and  4 fa!  show  the  time  outputs  of  the  two  lasers.  The  pulse 
width  of  che  mode-locked  output  was  about  l  nsec.  Figures  3(e)  and  4(e)  show  che  optical 
spectra  of  both  lasers  simultaneously.  The  resolution  of  the  optical  spectrum  analyser 
(-20  MHz)  was  not  sufficient  co  separate  out  che  two  spectra  for  typical  operation  where 
the  offset  between  che  two  spectra  was  about  10  MHz,  such  as  for  Figure  3(c).  However,  In 
Figure  4(e)  the  LO  offset,  Af,  was  set  (for  observational  purpoaes  only)  at  about  37  MHz 
so  chat  both  sets  of  modes  are  distinguishable,  i.e.,  each  alternative  peak  corresponds  to 
«  mode  of  che  seme  laser.  Each  laser  oscillates  in  about  20  longitudinal  modes  spread  over 
about  1.3  GHz.  Thus,  for  direct  observation  of  che  laser  signals  and  for  ranging,  che  de¬ 
tector  and  receiver  electronics  need  about  l.S  GHz  of  bandwidth.  Indeed,  the  lower  curve 
in  Figure  4(f)  shows  direct  ranging  results  using  fast  detectors  and  electronics  for  the 
two  point  scatterer  case. 

The  lower  trace  in  both  Figure  3(b)  and  Figure  4(b)  show  che  desired  filtered  time  sig¬ 
nal  from  che  target  after  heterodyning  with  che  LO  light.  Ftgurb  4(c)  shows  the  output 
after  a  rectification  stage  (compere  with  Figure  4(f)).  Notice  that  che  time  scale  la  2 
usec/dlv  in  Figure  4(c)  and  3  nsec/div  in  Figure  4(f).  The  comparison  shows  that  the  band¬ 
width  compression  f or  time  slowdown)  is  consistent  with  che  expected  ratio  of  F^/AF  *• 

300.  Notice  che  time  reversal  phenomenon  in  Figure  4(c).  Although  a  GHz  response  detector 
is  required  to  range  resolve  che  two  targets  in  Figure  4(f),  only  a  MHz  response  detector 
is  required  in  Figure  4(c). 

Figure  5 fa)  shows  part  of  che  detected  electrical  spectrum  for  a  single  point  target, 
similar  to  chat  shown  in  Figure  3(d)  but  on  an  ezpandad  scale.  Also,  the  mode-locked  dif¬ 
ference  frequency  AF  was  increased  to  338  kHz,  which  makes  che  individual  modes  more  evi¬ 
dent.  Figure  5(b)  shows  the  same  spectrum  but  displayed  using  an  acouscoopclc  (AO)  spec- 
crus  ahalvzer  setup  composed  of  a  1  mW  He-Ne  laser  beam  illuminating  a  Bragg  AO  cell  (-20 
MHz  bandwidth  centered  at  40  MHz),  a  lens  to  produce  the  FT,  and  a  1024  element  linear  de¬ 
tector  array  at  che  FT  plane.  The  detected  signal  from  the  DMML  setup  was  directly  fed 
into  che  AO  cell.  The  AO  generated  spectrum  for  che  two  point  target  case  (—16"  separa¬ 
tion)  is  shown  in  Figure  5(c).  The  expected  (-1  cycle)  sinusoidal  envelope  on  che  indi¬ 
vidual  cross-correlation  spectra  Is  clearly  evident.  Again  a  GHz  bendossa  phenomenon  is 
displayed  using  a  MHz  bandpass  dsvlce. 

Conclusions 

The  concept  of  reducing  the  bandwidth  requirements  on  che  detector  and  receiver  elec- 
t-onics  for  laser  ranging  applications  has  been  demonstrated  co  be  feasible  by  using  a 
mode-locked  transmitter  laser  and  a  mode-locked  local  oscillator  laser. 

This  work  was  suooorted  bv  the  U.S.  Army  Heseareh  Office. 
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Ftzura  1.  An  axamola  aae  of  spacera  for  a  3-ooda  lasar.  (a>  Tha  ear^ae  atonal  apacerum. 
<b>  Tha  local  oactllaeor  apacerum.  (cl  Tapfac  atznal  aueoeorralaeton  apacerum.  (d)  Local 
oactllaeor  aucocorralaeton  apacerum.  fa)  Tha  Xno  croaa-corralaelon  apacerum.  (f)  Tha 
X?q  croaa-corralaelon  apacerum.  (s'  Tha  coeal  oat ace ton  apaccrun.  (hi  Tha  oucpuc  fll- 
earad  apacerum. 


r’aure  2.  'Tha  axoarlraantal  s«cuo.  OSC-oaclllator ,  D-daeaceor,  BS-baamapUeeer,  L-lens, 
M-mttror,  and  PH-plnhola. 


-jl-  74  MH* 
FREQUENCY 


<c) 


J  j 


0  37  MHi 


FREQUENCY 

(d) 


!  FREQUENCY 

37  MHz  (•) 


-:a'jre  3.  Exoer f aental  r-sulc3  for  one  Do<nc  carzec.  (a'  The  :wo  He-Ne  node-locked 
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Figure  4.  Exoeritnencal  resit 
Ke-Ne  mode -locked  laser  out 
Cb'  Too  trace:  Che  differ 
Lower  trace.  the  received 
bandoass  filtering.  fc'  Th4 
cal  spectrum  of  the  detected 
3bout  27  MHz.  fe>  The'oot 
about  37  MHz.  fP  Uooer  tr 
for  direct  detection  with 
two  point  targets. 


Its  for  two  oolnc  targets  seoarated  by  14  Inchea.  (a)  The  two 
puts  acquired  via  fast  detectors  and  a  sampling  oscilloscope, 
once  signal  between  the  drive  signals  for  the  mode  lockers, 
signal  after  heterodyning  with  the  local  oscillator  light  and 
filtered  signal  after  rectification.  (d>  Parc  of  ths  electri- 
slgnal-  The  desired  soeccrum  to  be  filtered  is  centered  at 
ical  soectrum  of  both  lasers  when  the  LO  offset  frequency  is 
e:  the  transmitted  oulses.  Lower  trace:  this  received  signal 
ast  detectors  and  samollng  scooe  acquisition  for  Che  case  of 
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37  MHz 
FREQUENCY 

(o) 


I 

37  MHz 
FREQUENCY 

(b> 


37  MHz 
FREQUENCY 

(c) 


fMzure  Electrical  soectra  of  Che  detected  sicnal.  Horizontal  scalea  for  ail  Che 
Hzures  fa  2."  MHz'div  and  che  mode-lock  difference  frequency  la  5 38  kHz.  ,(al  and  (bl  show 
a  one  ooinc  carzac  resoonse  and  ^cl  show#  che  response  for  two  ooii.,  targets  with  a  15  inch 
seoaration.  (ai  shows  Che  resoonae  with  a  Tektronix  7L1 2  electronic  spectrum  analyzer. 
tb'  and  tc)  are  AO  cell  generated  aoectra.  Measured  resolution  of  the  AO  soeetrum  analyzer 
was  220  kHz. 
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APPENDIX  D 


TAPPED-FIBER  SPECTRUM  ANALYZER* 


C.  C.  Aleksoff  and  T.  M.  Stachelek 


Abstract 


An  optical  spectrum  analyzer  was  experimentally  demonstrated  that 
consisted  of  three  main  parts:  (I)  a  periodically  tapped  optical  fiber 
to  produce  a  sampled  one-dimensional  spatial  signal,  (2)  appropriate 
optics  to  produce  a  spatial  Fourier  transform,  and  (3)  a  holographic 
matched  filter  to  (In  effect)  phase  correct  all  the  optical  taps. 
Alternative  potential  configurations  and  signal  processing  operations 
are  also  presented. 


to  Optics  Communcatlons  for  publication. 


APPENDIX  E 


OPTICAL-HYBRID  BACKPROJECTION  PROCESSING* 


C.  C.  Aleksoff,  I.  J.  LaHale  and  A.  M.  Tal 


Abstract 


In  some  synthetic  aperture  systems*  the  data  gathered  from  a  single 
transmitted  pulse  samples  the  3-D  frequency  space  of  the  object.  A  3-D 
Image  can  be  formed  by  taking  a  3-D  FT  on  this  data.  However*  the  data 
are  often  In  a  form  that  allows  backprojectl on  techniques  to  be 
efficiently  used  to  generate  the  Image.  We  describe  an  optical  hybrid 
coherent  backprojectl on  processing  technique  that  uses  a  1-D  spatial 
light  modulator,  a  coherent  optical  processor*  a  2-D  detector  array,  and 
some  simple  computer  post-processing  to  produce  the  image. 


*Th1s  paper  was  presented  at  the  10th  International  Computing 
Conference,  "Unconventional  Imaging  and  Unconventional  Transformations" 
April  1983.  Published  in  the  Proc.  SPIE,  Vol  422,  1983*  pp.  89-95 
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OPTICAL— HY8R 10  8ACKPROJECTION  PROCESSING 


C.  C.  Aleksoff,  I.J.  Lartaie,  and  A.M.  Tal 


Environmental  Research  Institute  of  Michigan 
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Abstract 


In  some  synthetic  aperture  systems  the  data 
gathered  from  a  single  transmitted  pulse  samples 
the  3-0  frequency  space  of  the  object.  A  3-0 
image  can  be  formed  by  taking  a  3-0  FT  of  this 
data.  However,  the  data  are  often  in  a  form  that 
allows  backorojectlon  techniques  to  be  effi¬ 
ciently  used  to  generate  the  image.  We  describe 
an  optical  hybrid  coherent  backprojectlon  proces¬ 
sing  technique  that  uses  a  1-0  spatial  light 
modulator,  •  coherent  ootica!  processor,  a  2-0 
detector  array,  and  some  simple  computer  post¬ 
processing  to  produce  the  image. 

Introduction  -  Backprojectlon 


In  such  radar  and  laser  radar  systems  as 
spotlight  synthetic  aoerture  radar  (SSAR)  and 
those  that  image  rotating  objects,  the  construc¬ 
tion  of  the  3-0  image  U(x)  can  often  be  de¬ 
scribed.  as  will  be  shown  later,  by  the  3-0  FT 
(Fourier  transform)  f 1 —51 


i  /u If)  e"Z’,f”df 

where  U(x)  describes^the  (complex)  image  ampli¬ 
tude  ig  coordinates  x  •  (x,  y,  i)  of  the  image 
and  u(fl  describes  the  available  (camrlex,)  data 
in  coordinates  f  •  If„,  fy,  fz)  of  spatial 
frequency  space.  We  can  write 


a(f)  u0(f) 


where  u0ff)  represents  the  spatial  frequency 
descriotlon  of  the  ideal  image  U0(xl  over  all 
of  soace  and  a(»)  represents  the  3-0  aoerture 
over  which  the  data  was  gathered.  Often  the  data 
only  lies  on  some  3-0  surface,  as  will  be  shown 
Uter  for  some  specific  examoles.  In  any  case, 
the  actual  reconstructed  image  can  oe  described 
as  the  3-0  convolution  of  the  soread  function 
with  the  ideal  function  as  given  by 


U(x)  -  Ujj(x)  ***A(x) 


-9r3(«(n} 


is  the  3-0  screed  function  and  represents  a 
3-0  convolution.  We  note  that  an  active  area  of 
research  Is  the  extending  of  tr '  aperture  by 
filling  in  the  frequency  data.  a.g..  Interpola¬ 
tion  and/or  extrapolation,  but  our  only 
consideration  in  this  paper  Is  to  use  the  data 
directly  [6], 

Often  the  data  can  bn  described  as  being 
gathered  along  straight  lines  that  go  through 
the  origin,  which  we  will  call  polar  lines.  Let 
L(f,  ol  define  a  polar  line  (of  delta  width)  in 
the  direction  of  the  unit  vector  *  6.  Now  we 
note  that  the  3-0  FT  of  a  polar  line  Is  a  plane 
through  the  origin,  which  we  will  call  a  polar 
plane  [7,8].  Namely,  we  have 


P(x,  0)  .  t^{l(f,  b)J 


where  R(x,  b)  describes  a  polar  plane  with  delta 
support  in  the  direction  of  its  normal  b. 

Now  we  consider  the  case  where  the  data  In 
frequency  space  lies  along  polar  lines  as  can  be 
described  by  - 


u(0  -£u'(f)L(f,  ffl) 


u’(0  -  «#(f.)«*(f).  (8) 

Mere,  a* (f)  Is  the  part  of  the  aoerture  function 
describing  the  suooort  of  the  data  along  polar 
lines  orientated  in  the  direction,  fn,  T.e., 
the  tots!  aperture  description  is 


a(f)  -  a’ff)  Uf.  f„). 


I 


The  i'ir.a$e.  i.e..  the  FT  of  j(f)  as  given  by  £q. 
(71 ,  ca.i  be  written  as 

u( t)  .VtT’f x)  (i 


tT‘  '■  x  )  »  U'(x)***F(x,  x  )  (11) 

0  n  o  n 


-  f U«(5)  i(x  -  X  •  o)dq. 


ya  •  ya  *** 


*n  -  xxn*  fn  *  ffn  such  that 
,  xn  •  fn  ’  ' 


The  quantity  U'0  (•)  described  by  Eg.  (11)  is 
known  as  the  Radon  Transform  of  U * 3( • ) ,  and  is 
obviously  Just  the  integration  of  U'0  over 
planes,  i.e..  the  convolution  operation  gf  U'0 
with  P(-).  The  values  of  U'0  for  fixed  x„ 
is  described  as  orojections  of  U •  0( - >  '?]• 

Thus,  we  see  that  the  final  image  formed  with 
this  data  can  be  considered  as  the  sum  of  the 
Radon  transforms  of  the  ideal  image  convolved 
with  the  oartial  spread  function  A ■ ( - ) . 

Let  us  now  go  back  to  Eq.  (7)  and  derive  the 
operational  FT  form  that  we  seek.  Since  L ( * )  is 
a  line  delta  function  we  can  write 

u(f)  V«'(VL(?,  y.  (1*) 

n 

emphasizing  the  one  dimensional  nature  of  each 
line  of  data  by, using  fn  •  ffn,  where  it  is 
understood  that  fn  's  a  fixed  unit  vector, 
hence,  the  3-0  FT  of  u(7)  can  . then  be  written  as 

'J(x)  x)L(T,  xfl)  *  Pfx,  xp)  (15) 


xn  *  V  1 

The  quantity  vif - )  is  the  one  dimensional  FT  along 
the  nth  polar  line.  The  operation  described  by 
Eq.  (15)  is  known  as  backorojection  and  has  the 
slmole  geometrical  understanding  illustrated  in 
*  inure  1.  Each  1-0  FT  is  backorojected  along  a 
line  normal  to  <„  and  intersecting  the  desired 
'•construcflon  ooint  ~.  Thus,  we  can  write  very 
simoly  that 


Figure  1 .  The  Sackprojectlon  Technique 


.^Wfl(x  cos  »n) 


where  »n  is  the  angle  between  vector  to  the 
desired  reconstruction  point  1  and  the  polar  line 
direction  xn.  Thus,  it  is  only  necessary  to 
sum  backprojected  scaled  1-0  FTs.  He  will  demon¬ 
strate  this  type  of  orocesslng  optically  in  a 
later  section,  although  not  for  the  3-0  case  but 
rather  for  the  2-0  oolar  distributed  data. 
However,  the  conceot  is  viable  for  the  full  3-0 
case. 


Figure  2.  .'rage  Transform  Relationships 
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ffigure  2  snows  a  block  diagram  of  ways  to 
reconstruct  an  image.  In  coherent  systems  the. 
data  and  image  can  be  complex  whereas  in  inter¬ 
ferometric  or  tomograonic  systems  the  image  is 
nonnegative,  the  advantage  of  going  along  a 
backorojection  oath  rather  than  a  FT  oath  is  for 
comoutational  saving.  For  examole,  the  use  of 
the  FFT  in  digital  processing  requires  the  data 
to  be  on  a  oaral leoioed  grid,  which  imolies  that 
if  the  data  not  fall  on  the  grid  points  that 
interpolation  is  necessary.  The  interpolation 
comoutational  burden  can  often  exceed  the  FFT 
comoutational  requirements  by  many  orders  of 
magnitude  and  makes  the  FFT  based  technique 
unattractive. 

There  is  also  a  difference  as  to  where  the 
data  are  gathered.  In  tomography  the 
projections  are  the  basic  entity  gathered 
whereas  .in  syntnetic  aperture  systems  frequency 
space  data  are  usually  the  basic  entity  gathered. 


nthetic  aperture  Basics 


Consider  a  coherent  bi static  synthetic  aper¬ 
ture  imaging  system  that  illuminates  an  object 
or  scene  from  position  rg  and  receives  the  re¬ 
turn  at  position  r2.  See  Figure  3.  The  co¬ 
ordinate  system  is  assumed  fixed  with  respect  to 
the  object.  For  a  monochromatic  scalar  E-field 
'analytic  signal)  description,  received 
signal  from  a  differential  element  scatterer  at 
is  given  by  fwitMn  a  complex  constant 
multiplier) 


exp( ikR,)kR, 


exp{ ikR,)df 
i  0 

M31 


where  -fr01  is  the  scattering  coefficient, 
the  wave  number 


(19) 


u  the  temporal  frequency  and  R*  »[“,  -  FJ. 
Integrating  over  the  object  (all  the  scatterers) 
gives 


r  ik(R.rR-) 

*  j  —  d?0‘ 


In  the  far  field  we  can  let 

Ri  *  l?i  -  Fo|  *  r.i  '  V  V  f21) 

in  the  exponent  and  let  Ri  ■  rg  In  the  denom¬ 
inator  to  get 


2  ik(riT2)  r  ; 


'  i  2'  r 

— y 


dr0  (22) 


7  1  r  (rl  *  r2>  '  Z  (rl  *  r2>  (23) 


£d.  (22)  Is  a  basic  result  which  states  that  the 
received  signal  is  proportional  to  the  3-0  FT  of 
the  object  scattering  function  when  the  spatial 
frequency  is  properly  interpreted. 

We  would  like  to  invert  Eq.  (22)  to  find  a 
from  S.  In  the  cases  of  our  interest  the  in¬ 
version  is  helped  by  the  fact  the  frequency  v  is 
centered  about  some  carrier  frequency,  and  hence 
the  effect  of  k2  is  reasonably  constant  enough 
to  ignore.  Furthermore,  the  rate  of  change  of  v 
is  small  enough  to  allow  the  above  first  order 
(quasi-CW)  analysis.  Also,  rg  and  r2  do  not 
typically  vary  enough  to  have  to  be  taken  into 
account  other  than  in  the  exponent  where  the 
variation  of  k(rg  *  r?)  is  significant.  In 
fact,  to  account  for  this  factor  the  received 
signal  is  multiplied  by  a  function  proportional 
to  exp[-lk<H  *  r2)]  to  remove  its  effect 
and  truly  establish  what  is  known  as  a  MCP 
(motion  condensation  point)  in  space  fixed  with 
respect  to  the  object.  The  MCP  corresponds  to 
the  coordinate  center  of  our  frequency  space. 
Thus,  with  these  modifications  and  understand¬ 
ings,  we  can,  say  that  the  image  (estimate  of  the 
scatterers)  is  given  by 


(24) 

where  the  integral  is  evaluated  over  the  support 
of  the  available  data  points  in  frequency  space, 
i.e.,  over  the  3-0  synthetic  aperture.  Notice 
that  the  aperture  is  created  by  changes  in  the 
spatial  frequency  f  as  determined  by  the  param-  - 
eters  rg,  r2.  or 


4 


(a) 


butlon  in  frequency  space. 


Imaging  System  Examples 

Consider  a  SSAR  (Spotlight  Synthetic  Aperture 
Padart  where  the  aircraft  carrying  the  radar 
flies  along  a  straight  path'  and  keeps  the  radar 
oointed  to  the  same  are...  The  radar  transmits 
?*•  linear  chirp  waveforms  and  also  receives  the 
returns.  Thus,  Eq.  (241  holds  and  Eg.  (231 
becomes 

T  .  2rv/c  '  (251 

where  r  »  r-)  «  rj  for  this  case.  If  the  - 
radar  frequency  is  chirped  from  to 
then,  as  given  by  the  aPove  equation,  the  fre¬ 
quency.  vertor  direction  Is  stationary  but  the 
vector  magnitude  varies  from  fj  to  f?  where 
r",-  «  2?,v^/c.  Typ1cally,  the  chirp  pulses 
are  short  enouqh  for  r  not  to  change  signifi¬ 
cantly  due  to  the  aircraft  motion.  Thus,  we  see 
that  each  chiro  oroduces  a  section  of  a  solar 


f 

z 


(a) 


(b) 

Figure  5.  Imaging  a  Rotating  Object,  (a)  Data 
distribution  In  frequency  space  and  projection 
of  the  data  onto  the  f  ■  0  plane,  (b)  Re¬ 
sulting  image  slice  due  to  2-0  FT  processing 
projected  data  of  Fig.  5(a). 


line  of  data  In  frequency  soace  corresponding  to 
the  line  of  sight  direction  lo  the  motion  com¬ 
pensation  point,  as  indicate*  in  Figure  4.  All 
the  polar  lines  lie  in  a  plane,  called  the  slant 
plane,  since  the  radar  is  traveling  along  a 
straight  line.  The  general  1-0  problem  can  be 
considered  here  to  degenerate  into  a  2-0  pro¬ 
cessing  problem  since  all  the  data'  lies  in  a 
plane,  i.e..  the  aperture  is  in  the  slant  plane, 
and  hence  the  spread  function  is  Infinitely  long 
in  the  direction  perpendicular  to  the  slant 


best  a  2-0  projec- 
pto  the  slant  plane, 
(roly  a  2-0  FT  of 


plane.  Thus,  the  image  is  at| 
tion  of  the  3-0  ideal  image  o| 

The  required  processing  is  si 
the  polar  formatted  data. 

Let  us  consider  the  case  i>f  a  fixed  radar 
imaging  a  rotating  object.  Let  the  radar  trans-  ■ 
mlt  FM  linear  chirps  from  frebuancy  vi  to  v?. 

In  our  coordinate  system  fixed  to  the  object  the 
radar  appears  to  move  around  ihe  target  in  a  cir¬ 
cular  path.  Placing  the  motion  compensation 
point  on  the  object  rotation  4*1*  produces  odlar 


s 


lines  of  data  that  lie  on  the  surface  of  a  cone 
as  illustrated  in  Figure  5a.  The  finite  chirp 
truncates  the  cone  surface  between  radii  of  fj 
and  f?,  where  fj  »  ^(rj  *  r2>/c.  nie 
spread  function  for  a  truncated  cone  tends  to 
look  like  a  hollow  cone.  Typically,  the  3-0 
frequency  soace  data  can  be  projected  onto  the 
fz  »  0  olane,  producing  a  2-0  polar  formatted 
signal  wnich  can  be  processed  with  a  2-0  FT  to 
obtain  a  slice  fz  ■  0  plane)  of  the  object.  See 
Figures  5a  and  5b.  Scatterers  that  are  not  in 
the  slice  olane  will  produce  rings  instead  of 
points,  i.e.,  slices  of  the  spread  function  [11. 

Optical  Sackprojection  Processing 

In  the  previous  section  some  examples  were 
given  as  to  how  2-0  polar  formatted  inout  signals 
can  be  generated.  Figure  5  shows  the  classical 
technique  for  optically  processing  such  an  input. 
The  polar  lines  of  data  with  a  polar  offset 
spatial  frequency  carrier  are  written  onto  a  2-0 
spatial  modulator,  e.g.,  density  modulated  pt-cto- 
graphlc  film.  The  signal  is  then  coherently 
illuminated  and  the  FT  is  taken  optically  to 
produce  the  usual  conjugate  images.  One  of  the 
images  is  then  detected  for  further  analysis. 

In  its  simplest  terms  this  process  can  be 
described  by  writing  the  input  as 

vff)  •  8  *  2^Jun(f)j  cos  (2irf  •  h  ♦  #n) 

n  (26) 

where  3  is  the  bias,  h  the  polar  offset  vector, 
and 

Vf)  -  |u„w|  expf1dn(f)]  (27) 

are  just  the  individual  terms  of  the  summation 
given  on  the  right  hand  side  of  Eq..  (141.  Notice 
in  this  section  that  the  coordinates  and  vectors 
are  all  understood  to  be  2-0  rather  than  3-0  as 
in  the  earlier  sections.  After  the  FT  by  the 
lens  the  processor  output  amplitude  is  given  by 

V(x!  .  3«(x)  *^Un(x  -  R)  *  U*(x  *  h)  (281 


where  !Jn(*l  is  the  FT  of  uR(»)  and  *  reore- 
sents  conjugation.  We  can  identify  U„(»)  with 
the  terms  in  the  summations  in  £q,  (10)  or  Eq, 
'151.  That  is,  U„M  simoly  represents  the 
nth  backprojectlon  and  such  backprojectlons 
occur  naturally  in  this  system  due  to  the  input- 
formatting  and  the  FT  operation.  The  offset 
vector  h  f$  picked  to  seoarate  the  two  images 
away  from  the  zero  order.  Detection  of  one  of 
these  images  gives 
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Figure  6.  Classical  Technique  for  Optically  FT 
Processing  Polar  Formatted  Oata 


Memory 

□33X93 

Ref. 

source  Rotatable  ] 
/  detector 

Input  /  array — - 


!L  .  Of 

Backprojectlons1"' 

Figure  7.  The  Optical  8ackproject1on  Processor 

In  the  above  technique  all  the  Input  data, 
lines  are  2-0  FTed  simultaneously  in  parallel 
and  coherently  sunned  at  the  output  to  produce 
the  Image.  An  alternative  technique,  which  we 
call  the  backprojectlon  tefchnlque,  is  to  2-0  FT 
each  line  sequentially  (in  real  time)  and  then 
coherently  sum  them  in  a  memory  as  is  illustrated 
in  Figure  7.  The  detector  is  rotated  to  the 
proper  angle  for  each  input  line.  The  motivation 
for  this  alternative  technique  is  to  avoid  2-0 
spatial  modulators,  which  are  too  limited  in 
their  capabilities  for  our  use,  and  Instead  take 
advantage  of  the  newly  developing  technology  of 
spatial  detector  arrays  and  the  mature  technology 
of  1-0  modulators,  such  as  acousto-optic  cells. 

The  backprojectlon  technique  can  be  described 
by  letting  the  input  be  given  by 


v  (f)  -  8  *  Caff)  e’ 


U„(f)  cos  ( 2»f»h*#()) 


_  .  !  J 

where  Cjffle'v  is  the  reference  wave  (a  phase 
controllable  point  source  at  the  origin  of  the 
input  plane).  The  output  of  the  processor  (a 
2-0  FT)  is  given  by 


£d„{x  -  hi 


which  is  the  desired  result. 


*  Un(x  -  h)  ,*  U  (x •*  h)  (31) 


msmm® 


~he  detected  signal,  for  one  of  the  image  terms 
plus  the  reference  wave  is 

y*.  r>  -  jc  e'Y  *  y?  -  h > | 

,  I  -  12 


c  ‘  IV*  -  h1 


C  Un(x  -  h)  cos  (y  *  #n) 


Only  the  last  term  of  this  expression  is  desired, 
tt  can  6e  isolated  by  using  the  reference  wave 
Dhase.  Note  that 

I  'x,  01  -  tfx,  f)  .  2C  U„(x  -  h)  cos  in 

tv  “  n  n 


V*’  f1  -  V*-  -¥  *2C  Un(;-h)  s,".*n 

Thus,  we  can  obtain  the  quadrature  components  (I 
and  Q)  by  shifting  the  reference  wave  by  »/ 2  in¬ 
crements  for  each  line  input  and  then  subtracting 
the  correct  components.  The  subtraction  is  done 
digitally  and  the  results  are  stored  in  a  digital 
memory.  The  detector  (or  Input)  is  rotated  to 
the  next  position  and  the  procedure  is  repeated. 
The  running  sums  of  the  I  and  q  components  are 
keo:  and  the  image  at  any  time  can  be  calculated 
as  the  sum  of  the  squares  of  I  and  Q  to  give 


which  is  the  desired  result.  See  Eq.  (29). 

Experimental  Setup 

The  actual  experimental  implementation  dif¬ 
fered  from  the  conceptual  description  given 
above  and  tended  to  be  set  up  as  shown  in  Figure 
3.  The  experiment  was  performed  in  non-real 
time,  and  for  convenience,  the  1-0  data  lines 
were  generated  by  imaging  a  narrow  slit  onto  a 
2-0  ohotographic  recording  of  the  polar  formatted 
soatial  frequency  data.  The  input  data  and  image 
plane  were  held  stationary  relative  to  each 
other:  the  rotation  of  each  polar  slice  was 
achieved  instead  by  rotating  the  slit.  The  total 
processed  aperture. was  determined  by  the  number 
of  slices  coherently  summed  and  the  width  of  the 
slit,  both  of  which  were  allowed  to  vary.  The 
advantage  of  this  setup  was  that  by  fully  opening 
up  the  slit  a  classically  processed  coherent 
image  could  be  generated  and  directly  conpared 
to  the  image  processed  from  projections. 

The  image  plane  data  was  sampled  using  a  CIO  . 
array  camera.  The  required  subtractions  for  the 
in-onase  and  quadrature  (I  and  Q)  components  were 
performed  on  a  desk-too  comouttr  which  received 
the  digitized  data  from  the  CIO  array.  The  phase 
shifts  were  introduced  into  the  reference  wave 
by  monitoring  the  fringes  caused  by  the  inter¬ 
ference  of  the  reference  beam  and  the  undif- 
f 'acted  (OC)  portion  of  the  object  beam  with  a 
oair  of  detectors  whose  spacing  matched  that  of 
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Figure  8.  The  Experimental  Setup  to  Demonstrate 
Optical  Backprcjectton  Processing.  Notation: 

L  -  lens,  T  -  translator  Sff  -  stepper  motor, 

BS  -  beamsplitter,  K  -  mirror,  and  FT  -  fringe 
detector.  Solid  lines  are  electrical,  dashed  are 
optical,  and  dotted  are  mechanical  connectors. 


the  fringes.  The  detectors  were  operated  in  a 
difference  mode,  and  their  output  used  to  drive 
a  lock-in  mpllfier.  The  output  of  the  lock-in 
amplifier  was  used  to  control  a  piezoelectric 
mirror  in  the  reference  beam  path.  In  this  man¬ 
ner,  a  feedback  loop  was  created  which  served  to 
“freeze*  the  fringe  pattern  onto  the  detector 
pair,  i.e.,  to  keep  the  relative  phase  between 
the  object  and  reference  waves  steady.  By  moving 
the  detector  pair  laterally  across  the  fringe 
pattern,  this  “freezing*  effect  shifted  the 
reference  wave  phase  by  a  controllable  amount. 

The  same  feedback  loop  also  served  to  remove 
slowly-varying  fluctuations  in  the  relative  phase 
of  the  two  beans,  which  is  also  necessary  for 
proper  performance  of  the  experiment. 

The  data  acquisition  and  system  control 
functions  are  performed  electronically  by  both 
analog  and  digital  equipment.  The  overall  s'ystam 
is  controlled  by  a  Tektronix  4054  Graohics  System 
Computer.  Image  data  from  a  GE  TN2500  camera  in 
standard  interlaced-field  TV  format1 is  digi¬ 
tized  by  a  CVI  model  270  frame  grabber,  and  these 
are  supplied  *n  digital  form  to  the  4054  via  a 
dedicated  interface.  The  data  is  stored  by  the 
4054  In  an  auxiliary  memory  caoable  of  holding 
two  full  16-bit  grey  level  images  (2  x  256*  x 
16). 

The  system  reads  four  images  into  the 
graohics  conouter  for  each  position  of  the  aper¬ 
ture  slit,  with  relative  phases  of  0*.  180*,  90*. 
-90*  in  the  reference  beam,  and  perform*  the 
necessary  calculations  to  generate  tn#  I  and  Q 
images  containing  the  desired  information.  The 
aoerturt  is  then  rotated  by  a  specified  amount, 
and  the  procedure  repeated ,  Running  sums  of  the 
I  and  0  images  are  ktot  in  the  memory. 

The  data  collection  (i.e.,  image  formation 
processing)  did  not  run  in  real  time,  requiring 
aoproximately  45  seconds  per  subaperture.  An 
image  generated  from  64  projections  therefore 
requires  about  *5  minutes  of  run  time. 
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A  3asic  Result 

A  basic  result  is  presented  in  this  section 
for  a  single  point  target,  which  serves  to 
reoresent  the  "impulse  response"  of  the  back— 
projection  processing  algorithm. 

Figures  9a-9c  show  experimental  results  for 
the  reconstruction  of  the  coint  target,  tn  all 
cases,  the  length  of  the  slit  in  the  radial 
(range)  direction  is  13  mm.  Figure  9a  is  a  per- 
soective  plot  of  the  Imoulse  response  of  the 
slit,  which  in  this  case  has  a  width  of  0.2  mm. 
This  can  be  viewed  as  a  single  backprojection. 
Figure  9b  shows1  a  plot  of  the  image  obtained  by 
backorojection  processing  of  64  subaoertures. 

The  "comoression"  relative  to  Figure  9a  is  quite 
obvious.  The  effective  total  acerture  is  18  * 
12.3  me.  Figure  9c  is  a  plot  of  the  point  target 
image  produced  using  a  real  aperture  of  the  same 
dimensions. 

Conclusions 

It  was  shown  that  backprojection  processing 
can  be  a  viable  processing  technlgue  for  cases 
where  the  Fourier  freguency  space  of  an  object 
is  sampled  along  polar  lines.  Imaging  radar 
examoles  were  presented  that  were  characterized 
by  data  along  polar  lines.  The  concept  of  an 
optical-hybrid  backprojection  processor  was  suc¬ 
cessfully  demonstrated  experimentally.  The 
potential  for  using  a  1-0  spatial  light  modulator 
as  oart  of  a  real  time  system  makes1  the  technlgue 
very  attractive. 

This  work  was  supported  in  part  by  the  U.S. 
Army  Research  Office. 
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Figure  9.  Experimental'  Results  for  a  Point  Target, 
(a)  Single  backprojection.  (b)  Coherent  summation 
of  64  backprojections.  (c)  Classical  image  with  an 
aperture  equivalent  to  Fig.  9(b).  Horizontal 
scale  for  all  figures  is  the  same,  but  intensity 
(vertical  scale)  is  arbitral. 
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GRATING- BASED  INTERFEROMETRIC  PROCESSOR  FOR 
REAL-TIME  OPTICAL  FOURIER  TRANSFORMATION* 

Anthony  M.  Tal  and  Carl  C.  ATeksoff 


Abstract 


A  processing  approach  Is  Introduced  that  Is  capable  of  performing 
1-D  real-time  Fourier  transformations  on  the  Intensity  distribution  of 
an  Incoherent  optical  Input.  The  processing  approach  Is  based  on 
grating  Interferometers,  and  the  resulting  processors  are  simple  In 
structure  and  easily  Implemented.  Possible  processor  configurations 
together  with  experimental  results  demonstrating  the  operations  of  the 
system  are  presented.  Analyses  are  given  comparing  the  grating 
Interferometric  processor  to  the  Mlchelson  stellar  Interferometer  and 
the  classical  coherent  optical  processor. 
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Grating-based  interferometric  processor  for  real-time 
optical  Fourier  transformation 

Anthony  M.  TaJ  and  Cart  C.  Aleksoff 


A  processing  approach  ia  introduced  that  ia  capable  of  performing  1-D  real-time  Fourier  tranaformationa 
on  the  intenaity  diatribution  of  an  incoherent  optical  input  The  proceeding  approach  ia  baaed  on  grating 
interferometers,  and  the  resulting  processors  are  simple  in  structure  and  easily  implemented.  Poeaible  pro¬ 
cessor  configurations  together  with  experimental  results  demonstrating  the  operation  of  the  system  are  pre¬ 
sented.  Analyses  an  given  comparing  the  grating  interferometric  processor  to  the  Michelson  stellar  inter¬ 
ferometer  and  the  classical  coherent  optical  processor. 


L  Introduction 

The  Fourier  transformation  is  a  basic  operation  in 
many  signal  processing  applications,  and  it  has  been 
implemented  by  various  optical  techniques.  Coherent 
optical  processors1  can  perform  Fourier  transformation 
in  real  time,  but  they  have  several  familiar  limitations. 
Most  coherent  optical  processing  schemes  require  the 
.  use  of  an  incoherent-to-coherent  converter  which  re¬ 
duces  the  system  speed,  dynamic  range,  and  linearity. 
An  achromatic  Fourier  transform  system1-3  can  be  uti¬ 
lized  to  relax  the  temporal  coherence  requirement  and 
improve  the  system  SNR  by  reducing  coherence  noise. 
The  spatial  coherence  requirement,  however,  is  un¬ 
changed,  and  an  incoherent-to-  incoherent  converter  is 
still  needed. 

One-dimensional  processing  using  bulk  acoustical 
modulators  or  SAW  devices,4-5  on  the  other  hand,  are 
often  too  limited  in  their  space-bandwidth  product  by 
physical  constraints  such  as  transducer  frequency  re¬ 
sponse,  acoustical  velocity,  and  acoustical  attenuation. 
Moreover,  processing  with  an  acoustical  modulator  is 
limited  to  serial  inputs.  To  be  able  to  process  parallel 
spatial  inputs  such  as  images,  a  spatial-to-temporal 
conversion  is  needed. 

Incoherent  optical  processing  techniques  such  as 
OTF  synthesis8-7  permit  the  use  of  incoherent  optical 
inputs,  but  their  applications  are  limited  to  a  rather 
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restrictive  class  of  operation.  For  example,  complex 
Fourier  transformation  has  not  been  demonstrated. 
Cosinusoidal  transformation  was  demonstrated  by 
Mertz  and  Rogers8  using  a  shadow  casting  technique. 
Such  a  technique  operates  in  the  geometric  optics  re¬ 
gime,  and  the  achievable  space-bandwidth  is  conse¬ 
quently  sraalL  The  optical  vector-matrix  multiplier9-10 
can  be  configured  to  perform  Fourier  transformation; 
unlike  other  more  conventional  optical  processors,  the 
vector-matrix  multiplier  operates  on  discrete  sampled 
signals.  The  difficulties  in  constructing  and  aligning 
a  large  2-D  mask  and  the  limits  on  the  densities  of  these 
masks  restrict  the  space-bandwidth  product  that  can 
be  achieved  in  practice.  The  relative  complexity  of  the 
system  also  diminished  its  attractiveness. 

In  this  paper,  we  introduce  an  optical  processor  which 
belongs  to  a  class  of  processor  which  we  call  interfero¬ 
metric  processors.  This  processor  operates  on  inco¬ 
herent  optical  inputs,  but  its  performance  is  closely 
related  to  coherent  optical  processing  systems.  An 
achromatic  grating  interferometer  is  used  in  our  im¬ 
plementation  of  the  processing  approach.  It  is  capable 
of  producing  and  displaying  in  real  time  the  real  and 
imaginary  parts  of  the  Fourier  transform  of  a  1-D  in¬ 
coherent  optical  input.  Preliminary  experimental  re¬ 
sults  are  presented  to  verify  the  theory  and  to  demon¬ 
strate  the  processing  approach. 

H.  Complex  Degree  of  Spatial  Coherence 

Consider  the  system  geometry  illustrated  in  Fig.  1. 
The  radiation  emitted  by  a  space-limited  object  Held 
at  the  («, /3)  plane  is  detected  at  far-field  locations 
(xiO'i)  and  (x-j.jyi;)-  If  the  radiation  is  quasi-mono- 
chromatic  and  spatially  incoherent,  it  can  be  shown  that 
the  complex  degree  of  spatial  coherence  (CDSC)  is  equal 
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Fig.  1.  System  g tom* try  relating  the  intensity  distribution  of  an 
incoherent  object  Held  to  the  complex  degree  of  spatial  coherence  at 
a  detection  plane  in  the  far  field. 


C(u.u)  ■  K  KaJ})  exp  |  (ua  +  u£)jdad0,  (1) 

where 

um(x  t~xi)/R,  c- (>i -yj)/R,  1 
and  K  is  a  constant. 

We  see  that  there  is  a  Fourier  transform  relationship 
between  the  input  intensity  distribution  and  the  CDSC 
at  far  Held.  This  relationship  is  sometimes  known  as 
the  van  Cittert-Zernike  theorem.  Thus,  to  obtain  the 


frequencies  are  sampled  by  physically  varying 
separation  of  the  detecting  apertures.  Such  a  sequen¬ 
tial  measurement  scheme  is  too  restrictive  for  use  in  a 
real-time  processing  system.  The  requirements  that 
the  input  be  quasi-monochromatic  also  makes  the  sys¬ 
tem  very  inefficient  for  inputs  that  are  naturally  poly¬ 
chromatic. 

Roddier  et  al. 12  introduced  a  rotation  shearing  in¬ 
terferometer  for  astronomical  imaging.  Basically  the 
input  wave  front  is  separated  into  two  parts.  One  part 
is  rotated  180°  with  a  system  of  roof  prisms  or  dove 
prisms  and  then  recombined.  At  the  center  of  the  field 
which  is  the  center  of  rotation  there  is  no  shear.  Mov¬ 
ing  away  from  the  center,  increasing  amounts  of  shear 
are  introduced.  Thus,  unlike  the  Michelson  stellar 
interferometer,  all  the  spatial  frequencies  can  be  dis¬ 
played  at  once.  This  technique  is  attractive  for  many 
applications,  and  it  was  suggested  as  an  alternative  to 
stellar  speckle  interferometry.  However,  similar  to  the 
Michelson  stellar  interferometer,  it  requires  quasi¬ 
monochromaticity  and  is,  therefore,  rather  inefficient 
in  the  processing  of  inputs  that  are  spectrally  wide  band. 
George  and  Wang13  combined  the  rotation  shearing 
interferometer  with  an  achromatic  transform  optical 
system  and  demonstrated  polychromatic  cosine  trans¬ 
form.  Their  system,  however,  requires  a  system  of 
prisms  to  perform  the  rotation  and  a  multielement  op¬ 
tical  system  for  the  achromatic  transformation. 
Moreover,  switching  from  cosine  to  sine  transformation 
cannot  be  easily  achieved  due  to  the  difficulty  in  real¬ 
izing  achromatic  t/2  phase  shifts  using  reflective  and 
refractive  optics. 

In  this  paper,  we  introduce  a  grating-based  achro¬ 
matic  optical  processor.  It  is  capable  of  producing,  in 
real  time,  cosine  and  sine  transformations  of  inputs  that 
are  spatially  and  temporally  incoherent.  One  major 
advantage  offered  by  grating-based  systems  is  sim¬ 
plicity.  The  grating  interferometric  processor  tends  to 
require  much  less  hardware  than  conventional  systems 
making  it  easier  to  fabricate  and  set  up.  The  simplicity 
of  the  grating  system  also  makes  the  processor  more 
compact  and  vibration  resistant. 


Fourier  transform  of  the  intensity  distribution  of  an 
incoherent  input,  one  merely  has  to  devise  a  means  to 
measure  the  CDSC  at  far  field;  This  can  be  achieved 
interferometrically  by  measuring  the  complex  visibility 
of  the  interference  fringes  formal  between  the  fields  at 
(xiO'i)  and  (x20'2)-  The  fringe  visibility  function  is 
equal  to 

/UlO'il  + 

•  (21 

For  a  homogeneous  incoherent  input,  /(x iO'i)  * 
/ and  the  CDSC  is  space  invariant  in  the  far  field. 
Thus  we  may  write 

V'(«a,l  C(uu->  Cll 

With  the  Michelson  stellar  interferometer  and  other 
similar  interferometric  imaging  systems,  the  spatial 


W.  Achromatic  Grating  bitartaromatar 

The  achromatic  grating  interferometer  has  been 
studied  extensively  and  applied  to  many  applications. 
Chang  and  Leith14*16  applied  the  interferometer  to  in¬ 
terferometric  imaging  and  nondestructive  testing. 
Leith  and  Roth17-18  studied  the  noise  performance  of  the 
interferometer  and  the  synthesis  of  the  convolution 
integral  for  linear  filtering.  Collins19  demonstrated  the 
construction  of  holograms  and  matched  filters  with 
incoherent  light  Leith  and  Swanson20  utilized  the 
interferometer  for  the  fabrication  of  noise-free  dif¬ 
fractive  optical  elements.  Tai  and  Aieksoff21  employed 
the  grating  interferometer  arrangement  in  an  imaging 
system  for  improved  imaging  through  scattering 
media. 

Leith  and  Chang,10  in  particular,  have  recognized  the 
ability  of  the  grating  interferometer  to  perform  Fourier 
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transformation*.  Their  interest,  however,  was  pri¬ 
marily  in  imaging  where  the  image  of  an  object  was  re¬ 
constructed  by  performing  an  optical  transform  with 
a  coherent  optical  processor.  To  be  able  to  utilize  the 
interferometer  as  a  processor,  an  appropriate  decoding 
technique  must  be  developed  to  produce  a  complex 
output  that  can  be  efficiently  read  out  or  interfaced  to 
a  computer. 

The  simplest  achromatic  grating  interferometer  that 
can  be  realized  is  composed  of  a  single  grating.  Con¬ 
sider  the  optical  arrangement  depicted  in  Fig.  2(a).  For 
a  single  point  source  located  at  at »  ao,  the  field  ampli¬ 
tude  in  the  volume  where  the  ±1  diffracted  orders 
overlap  is  equal  to 

E(at)ixj)  ■  Ci/5(a«j 

•  |«sp(i2*r(Pe  -F)x  +  Huy/  1  -  -  FP) 

+  ttp(i2r(/*o  +  F)x  f  ifciy'  1  -  +  TPJ),  (4) 

where  k  ■  2r/X,S(ao)  is  the  intensity  of  the  point 
source,  C  is  a  constant,  and  Pq  “  sin(tan-lao//)/X,  /  is 
the  food  length  of  the  lens,  F  is  the  spatial  frequency  of 
the  grating,  and  X  is  the  wavelength  of  the  source.  The 
corresponding  light  intensity  distribution  in  the  fringe 
volume  can  be  written  as 

Kwa)  -  |C|  |l+co*  jtvf*  +  ” 


■  (vr-w« +r>*  - 


Fif.  2.  Sinck-sratioc  •chromatic  imarfaromater:  (a)  bsiic  geom¬ 
etry;  (b)  usins  a  baam  ipiittar  and  a  taiaacopk  imafins  ayawm  to 
daeoda  tbs  piWMaor  output:  (c)  imafinc  tha  ftinga  voiume  to  a  piana 
away  how  tha  gratins  to  aocaaa  tha  aaro  ipatial  frequency. 


a  spatially  incoherent  line  source  with  intensity  distri¬ 
bution  f(a )  extending  from  a  «  0  to  a  »  1,  the  output 
intensity  distribution  within  a  constant  factor  is 


Hxjt)  -  /(a)  Jl  +  coa4xF  |x  -  ^j|  da. 

The  intensity  along  the  optical  axis  at  x 

/(*)  ■  /<«*><!  +  coa matda 

m  So  fia)da*  f',M 


17) 


Ois 


Icoaw ada. 


18) 


(5) 


Using  a  first-order  approximation,  the  intensity  dis¬ 
tribution  can  be  expressed  ss 

-  |C|*S(o#)  |l  +  eoa  [srf  j*  -  y||  ■  (0) 

Now,  if  instead  of  having  a  single  point  source,  we  have 
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where  w  ■  ArFz/f. 

The  first  term  is  simply  a  constant  bias,  and  the  sec¬ 
ond  term  we  recognize  as  the  cosine  transform  integral 
Thus  by  reading  out  the  intensity  distribution  along  the 
optical  axis,  we  obtain  the  bipolar  cosine  transform  of 
the  1-D  input  intensity  distribution  plus  a  constant  bias. 
We  emphasize  that  the  expression  in  Eq.  (8)  is  not  a 
function  of  wavelength.  The  transformation  is  per¬ 
formed  achromatically. 

The  single-grating  interferometer  provides  an  ex¬ 
tremely  simple  means  of  achieving  optical  transfor¬ 
mation  achromatically.  However,  the  information  is 
displayed  as  an  intensity  distribution  along  the  optical 
axis,  which  is  rather  inaccessible.  In  particular,  the  zero 
spatial  frequency  is  located  at  z  »  0,  right  at  the  surface 
of  the  grating.  One  way  to  read  out  the  intensity  dis¬ 
tribution  along  the  z  axis  is  to  place  a  very  thin  beam 
splitter  along  the  z-y  plane  at  x  ■  0  as  shown  in  Fig. 
2(b).  The  intensity  distribution  can  then  be  imaged  to 
a  convenient  output  plane  to  be  measured.  Unfortu¬ 
nately,  it  is  not  straightforward  to  place  such  a  beam 
splitter  light  up  against  the  grating.  It  nevertheless  can 
be  done  using  special  optics.  One  may  image  the  fringe 
volume  directly  with  a  telescopic  imaging  system  as 
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Fig.  3.  Multiple-grating  achromatic  interferometer. 


shown  in  Pig.  2(c)  to  make  the  fringe  volume  more  ac¬ 
cessible.  However,  the  diffraction  angles  would  be  se¬ 
verely  limited  by  the  numerical  apertures  of  the  lenses. 
A  better  approach  is  to  utilize  a  three-grating  arrange¬ 
ment  similar  to  that  of  Leith  and  Chang14  as  shown  in 
Fig.  3.  The  first  grating  G  t  has  a  spatial  frequency  of 
F,  the  second  and  third  gratings  are  copianar,  and  both 
have  a  spatial  frequency  of  2 F.  The  input  light  field  is 
split  into  two  first  orders  and  the  first  grating  and  then 
recombined  by  the  second  and  third  gratings,  G  2  and 
G 3.  A  new  fringe  volume  is  formed  about  2  ■  2D,  ac¬ 
cessibly  away  from  any  grating.  Moreover,  spurious 
diffraction  orders  can  now  be  easily  removed  to  elimi¬ 
nate  the  spurious  fringes  and  added  bias  that  may 
otherwise  be  introduced.  With  a  single  point  source 
located  at  a  -  a0,  the  light  field  at  the  new  fringe  vol¬ 
ume  is  equal  to 

Blew r*y>  *  Cy/Siao)  ■  |eip(i2T(f>a  -  F)x 
+  i*Ov  1  -  A^  A)  +  F>2 
+  ik(t  +  D>vl  -  A2(Aj  -  F)-| 

+  txa[i2x{Po  +  F)x  +  ikD\  1  -  -  /)*■ 

+  ik(x  +  0>v  I  -  +  F)-’||.  (9) 

The  corresponding  light  intensity  distribution  can  be 
written  as 

Hxj)  m  |C|*S(a)(l  +  cctHwFx 

+  *i(vl  -  .\J(fo  +  75* -  V  1  -  A -</>„-  uo> 

which  can  be  approximated  by 

Hxj)  »  (C^Sfal  jl  +  coa  “  ~j||  ’  1111 

For  a  1-D  signal  along  the  a  axis  ranging  from  0  to  /, 
the  output  intensity  distribution  becomes 

/<  xjt  *  fiat  |l  +  cosJteF  jr  -  — jJJ  da.-  (12) 

At  x  *  0.  once  again  we  have  the  cosine  transform  re¬ 
lationship 

/It)  ■  /lift|  1  +  cixiwi' 'Id it.  1 13) 

To  perform  a  Fourier  transformation,  both  real  and 
imaginary  output  data  are  required.  By  simply  trans¬ 
lating  the  grating  (#>  along  the  .V.  direction  by  one- 
fourth  of  a  period  li.e..  t/2  phase  shiui.  the  light  in¬ 
tensity  distribution  along  the  ;  axis  becomes 


/(*)  *  f(au(t  +-sin(«ia)|<<a,  (14) 

which  is,  of  course,  the  sine  transform  of  the  input  f(a). 
Thus  both  the  real  and  imaginary  parts  of  the  complex 
transform  output  can  be  obtained.  We  may  note  that 
the  shifting  of  the  grating  produces  a  t/2  phase  shift  for 
all  wavelengths.  The  achromaticity  of  the  system  is 
maintained.  This  is  in  contrast  to  other  interferometric 
systems  using  reflective  or  refracting  optics. 

Once  again,  the  intensity  distribution  along  the  z  axis 
can  be  read  out  by  placing  a  thin  beam  splitter  along  the 
x-z  plane  as  shown  in  Fig.  4.  The  in tensity  distribution 
along  the  beam  splitter  is  mapped  onto  a  convenient 
output  plane  outside  the  fringe  volume  using  a  one- 
to-one  telescopic  imaging  system. 

This  simple  implementation  causes  several  problems. 
The  beam  splitter  has  to  be  very  thin;  a  pellicle  beam 
splitter  might,  therefore,  be  used.  Pellicle  beam 
splitters  unfortunately  have  a  tendency  to  vibrate. 
Although  this  is  not  a  fundamental  problem  and  can  be 
solved  using  special  optics,  we  decided  to  use  an  alter¬ 
nate  arrangement  for  the  demonstration. 

The  light  distribution  in  thenutput  fringe  volume  is 
equal  to 

Uxjrxl*  /(a)  jl  +  coa  (x -^jjjda. 

If  we  place  a  fourth  grating  with  spatial  frequency  F  at 
2  **  2o  as  shown  in  Fig.  5,  the  grating  will  demodulate  the 
light  back  to  zero  spatial  frequency,  and  the  output  light 
intensity  becomes 

Utal  •  J*  Ha)  jl  +  coa  j— jj  da ,  ( 15) 

which  in  effect  gives  the  cosine  transform  of  {(a)  sam¬ 
pled  at  spatial  frequency  4xFzo/f.  By  translating  the 


Fif.  4.  Mapping  the  transform  output  in  tho  fringt  volume  to  a 
convenient  output  plane  using  a  beam  splitter. 


Fig.  V  Simpling  -lutiai  treuuemie*  with  a  translating  fourth 
grating. 
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which  can  be  written  as 


/(«')  ■  f  /(oWa  +  /(a)  codw'aMa,  (17) 


where 


«'  «  (4»fy  tantl //, 

and  0  is  the  tilt  angle  of  the  fourth  grating. 

Once  again,  a  biased  output  of  the  cosine  transform 
is  obtained.  However,  the  spectrum  is  now  displayed 
along  the  w'  direction,  which  is  orthogonal  to  the  a  axis. 
The  cosine  spectrum,  therefore,  lies  along  a  direction 
that  is  perpendicular  to  the  1-0  input. 

If  we  translate  the  fourth  grating  by  one-eighth  of  a 
period  (i.e.,  ±r/4  phase  shift),  we  obtain 

/•i  /*< 

/(«*')  ■  |  f{a)da  ♦  /(<*>  uniw'ak/ar  (19) 

or  the  sine  transformation  of  /(a). 

To  perform  complex  Fourier  transformations,  we 
need  both  the  real  and  imaginary  parts.  The  cosine  and 
sine  transforms  can  be  displayed  simultaneously  by 
using  a  split  grating  as  the  fourth  grating  at  the  output 
as  shown  in  Fig.  7.  The  grating  lines  at :  he  two  halves 
of  the  grating  are  phase  shifted  by  *74.  This  produces 
a  cosine  transform  at  one  tide  of  the  output  and  the  sine 
transform  at  the  other. 

We  note  that  the  tilted  grating  decoding  technique 
may  also  be  applied  to  the  single-grating  interferometer 
by  placing  the  tilted  grating  against  the  first  grating 
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transvening  the  fringe  volume.  However,  the  presence 
of  spurious  diffraction  orders  may  substantially  degrade 
the  transform  output.  In  addition,  transform  infor¬ 
mation  near  dc  may  not  be  displayed  due  to  blockage 
by  the  edge  of  the  tilted  grating. 

When  only  cosine  transforms  are  performed,  the 
input  must  be  either  single  sided  (i.e.,  0  <  a  <  l)  or  an 
even  function  [i.e.,  /(or)  *  /( — «))  to  avoid  ambiguity. 
With  a  dual-channel  system,  however,  the  input  is  no 
longer  limited  to  being  single  sided  or  even.  The  pro¬ 
cessor  is  actually  performing  the  decomposition 


Fig.  6.  Decoding  the  immform  output  with  a  tiitad  fourth  grating. 


fourth  grating  along  the  z  axis  over  the  entire  fringe  box, 
the  cosine  transform  output  can  be  scanned  out. 
However,  such  a  translation  technique  takes  time,  and 
our  goal  is  real-time  operation. 

We  note  that  Sq.  (15)  ia  not  a  function  of  z  and  y. 
The  sampled  cosine  output  can  be  obtained  by  mea¬ 
suring  the  intensity  at  any  point  of  the  output  plane.  If 
we  tilt  the  fourth  grating  toward  the  x-z  plane  so  that 
it  traverses  the  fringe  volume  as  illustrated  in  Fig.  6, 
different  parts  of  the  grating  along  the  y  direction  will 
be  sampling  the  intensity  at  different  positions  along 
the  z  direction.  The  output  intensity  distribution  be- 


J*  /(<*)e*p(— iua)da  < 


X> 


eca<u>a)rfa 


-  *  J*  /(a)  sin(w aWa 


where  f(a)  is  an  arbitrary  real  and  non-negative  func¬ 
tion.  which  in  general  is  a  combination  of  even  and  odd 
parts.  That  is,  the  decomposition  ia  not  in  terms  of  the 
cosine  transform  of  the  even  part  and  the  sine  transform 
of  the  odd  part  but  rather  in  terms  of  7,  and  7/,  which 
are  the  real  and  imaginary  parts,  respectively,  of  the 
Fourier  transform.  In  this  operating  mode,  the  input 
can  extend  from  —/to/  with  unique  information  on  both 
sides,  and  the  space-bandwidth  product  of  the  proces¬ 
sor  is  in  effect  doubled.  This  doubling  of  the  space- 
bandwidth  product  can  be  considered  to  occur  because 
the  system  acts  as  a  dual-channel  processor  simulta¬ 
neously  performing  the  real  and  imaginary  parts  of  the 
Fourier  transformation.  A  more  detailed  discussion  on 
the  space-bandwidth  product  of  the  grating  interfero¬ 
metric  processor  is  given  in  a  later  section. 

IV.  Experimental  Results 

An  interferometric  optical  processor  was  set  up  using 
the  configuration  illustrated  in  Fig.  6.  As  the  input,  a 
ground  glass  diffuser  was  backilluminated  by  an  unfil¬ 
tered  xenon  arc  lamp  to  produce  an  incoherent  light 
source,  which  in  turn  was  used  to  illuminate  a  slit  and 
transparency.  To  test  if  the  system  can  produce  the 
correct  Fourier  transform  output,  a  transparency  with 
an  intensity  transmittance  of  1  +  cos(2*-ga)  was  used. 
Thuc  the  input  is  given  by  {(a)  -  rect((a  -  «o)f«|[l  + 
cos(2rga)|.  The  cosine/sine  transform  outputs  were 
read  out  using  a  1024-element  linear  detector  array.  In 
Fig.  8(a),  we  show  the  cosine  transform  output  obtained 
with  the  interferometric  processor.  For  comparison, 
we  computed  the  Fourier  transform  of  a  function: 


L^~ 


"&// 


Fif.  7.  Si  multanaoua  cotin*  and  tint  tranafnrmatiom. 


obtained  and  digitally  computed  results.  The  small 
phase  error  might  have  been  caused  by  irregularity  in 
m  Opoi  com*  the  local  frequency  of  the  gratings. 

The  fourth  grating  was  then  translated  by  one-eighth 
of  a  period  to  produce  a  jt/4  phase  shift  on  the  recom¬ 
bined  beams.  The  optically  obtained  sine  transform 
output  is  shown  in  Fig.  8(c).  For  comparison,  the  cor¬ 
responding  computed  imaginary  part  of  the  Fourier 
transform  is  shown  in  Fig.  8(d).  Once  again,  the  optical 
output  matches  the  digitally  computed  result. 

Mo^atc^a  V.  Comparison  to  a  Mchataon  Stellar  Interferometer 

The  operation  of  this  system  can  be  explained  more 
intuitively  by  going  back  to  the  Michelson  stellar  in¬ 
terferometer.  The  Michelson  stellar  interferometer  in 
its  basic  form  is  simply  two  apertures  separated  by  a 
distance  d  as  shown  in  Fig.  9(a).  The  CDSC  is  obtained 
by  measuring  the  complex  visibility  of  the  fringes 
formed  for  different  separations  of  the  apertures.  We 
m  c*iai  sa  fr0m  Fig.  9(b)  that  we  have  a  similar  situation  with 

the  grating  interferometer.  The  interfering  light  fields 
are  sheared  by  an  amount  Ax  mz*  tan(sin-1(X/)),  where 
F  is  the  spatial  frequency  of  the  grating.  Thus,  similar 
to  the  Michelson  stellar  interferometer,  two  parts  of  a 
wave  front  separated  by  a  transverse  distance  Ax  are 
made  to  interfere,  and  the  resulting  fringe  visibility 
along  the  z  axis,  therefore,  follows  a  Fourier  transform 
relationship  with  the  source  intensity  distribution. 
Unlike  the  Michelson  interferometer  where  the  visi¬ 
bility  for  each  value  of  Ax  must  be  measured  separately, 
the  grating  interferometer  presents  the  interference 
fringes  for  the  different  shear  distances  simultaneously 
at  different  positions  of  x.  This  interferometric  optical 

Pit-  a  Optically  CMMratedcoeine  and  line  tranat'unns  vs  the  digitally 
computed  versions. 


/'(a)  m  Wia\  rect  |1  ♦  ciwi^xeit  ♦  *J|  tipitd).  <  19) 

where  rect((«  —  «o)/£.|  “  1  for  <*«  -  L,  2  <  <t  £  <to  +  L! 2, 

W(a)  is  an  apodizing  factor  added  to  account  for  the 
nonuniform  illumination  of  the  input.  «.>  is  the  offset 
of  the  input  from  the  axis,  g  is  tne  spatial  frequency,  <fi 
is  the  phase  of  the  input  cosine  function,  and  0  is  an 
arbitrary  phase. 

The  choice  of  the  grating  position  that  is  designated 
as  cosine  of  sins  transform  is  arbitrary  and  unimportant 
as  long  as  an  exact  phase  shift  is  introduced  between 
the  sine  and  cosine  transforms  and  a  consistency  is 
maintained.  The  difference  in  the  outputs  caused  by 
different  initial  choices  of  the  grating  position  is  simply 
a  constant  phase.  (A  similar  situation  arise*  in  coherent 
optical  processing  where  a  coherent  reierence  beam  is 
introduced  to  measure  the  complex  output  of  the  optical 
processor.)  In  our  experiment,  the  constant  phase  term 
exptid)  was  obtained  by  varying  it  in  the  computation 
until  the  real  part  of  the  zero  order  of  the  digital  Fourier 
transform  matched  that  of  the  optical  output.  In  Fig. 

81b),  we  shuw  the  computed  real  part  o|  the  Fourier 
transform.  Except  for  a  small  phase  diiterence.  there 
is  a  very  good  correspondence  .between  the  optically 
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Fi*.  10.  Shifting  of  fring*  volume  for  off-axi*  input  point 


processor  is,  therefore,  basically  an  interferometer  that 
is  capable  of  measuring  in  parallel  the  1-D  CDSC  of  a 
field  emanated  from  a  homogeneous  incoherent  line 
source  at  far  Held  as  created  by  the  collimating  lens. 

The  spatial  frequency  as  measured  by  a  Michelson 
interferometer  is  given  by  u  *  2xA x/\f,  which  is 
wavelength  dependent.  A  Michelson  interferometer 
can,  therefore,  operate  only  with  quasi  -  monochromatic 
radiation.  To  make  the  system  achromatic,  a  variable 
amount  of  shear  must  be  introduced  so  that  Ax(A)/A  * 
constant.  The  amount  of  shear  introduced  by  the 
grating  interferometer  varies  with  wavelength,  since  to 
the  first  order.  Ax  ■  2FAz,  as  illustrated  in  Pig.  9(c). 
The  spatial  frequency  is  thus  u  »  4  rcFzlf,  which  is 
wavelength  independent  and  hence  achromatic. 

VL  Comparison  to  Coherent  Optical  Processing 

Although  the  interferometric  processor  operates  on 
incoherent  inputs,  its  performance  characteristics  are 
very  similar  to  a  coherent  optical  processor.  For  a  co¬ 
herent  system,  the  transform  relationship  can  be  written 
as 


F<P>  •  J* . /<“• **p (‘ ~\7^j  da • 

and  for  the  interferometric  optical  processor, 

F*U)  -  J*  f(a)  co*  da 

and 


(20) 

(21) 


F/(*l  ■  J*  /(a)  sin  |4*y*°j  da.  (22) 

The  difference  between  the  two  processors  is  in  the 
scaling  factors  2x/A/  and  4 xF/f.  The  scaling  factor  for 
a  coherent  optical  processor  is  inversely  proportional 
to  wavelength,  while  the  scaling  factor  for  the  interfer¬ 
ometric  processor  is  independent  of  wavelength.  Hence 
the  interferometric  processor  is  equivalent  to  an  ach¬ 
romatic  Fourier  transform  system.--  ' 

,  For  systems  that  use  AO  modulators  or  discrete 
sources  (e.g.,  the  vector-matrix  multiplier)  as  input 
devices,  the  space-bandwidth  product  is  typically  input 
device  limited  and  can  be  easily  defined.  For  a  coherent 
optical  processor  and  an  interferometric  optical  pro¬ 
cessor,  on  the  other  hand,  the  space-bandwidth  product 
is  usually  processor-limited  and  is  dependent  on  many 
factors  such  as  the  exact  optical  configuration,  the  //No. 
of  the  transform  or  collimating  lens,  the  aperture  of  the 
lens,  and  the  performance  of  the  lens,  in  the  following. 


we  compute  the  space-bandwidth  product  (SBWP)  of 
the  single-grating  and  toe  multiple-grating  interfer¬ 
ometers.  The  approximations  and  assumptions  used 
in  the  computation  are  listed: 

(1)  Hie  angular  bandwidth  of  the  gratings  in  the 
processor  are  very  wide  (e.g.,  surface  relief  gratings). 

(2)  The  spurious  diffraction  orders  ars  not  a  problem, 
and  no  stops  are  used  in  the  system.  (This  assumption 
is  always  valid  for  systems  that  utilize  beam  splitters  to 
demodulate  the  output.  It  is  also  valid  for  the  multi¬ 
ple-grating  system  using  a  tilted  fourth  grating  for 
output  demodulated  if  the  gratings  are  properly  fabri¬ 
cated.  The  effects  of  spurious  diffraction  orders  will 
be  discussed  in  a  later  section.) 

(3)  The  gratings  are  of  sufficient  size  that  the  entire 
field  over  the  full  angular  bandwidth  is  diffracted  by  the 
gratings. 

(4)  The  field  distribution  along  the  optical  z  axis  is 
considered  as  the  output  plane.  [In  actual  system  op¬ 
eration,  it  is  necessary  to  transfer  this  output  plane  to 
a  more  accessible  plane  using  a  beam  splitter  or  a  grat¬ 
ing.  This  analysis  is  directly  applicable  to  a  system  that 
utilizes  a  beam  splitter  to  demodulate  the  output  (Fig. 
4).  With  the  use  of  a  tilted  fourth  grating,  however,  the 
performance  characteristic  may  be  degraded  in  the 
decoding  process,  and  the  system  SBWP  may  be  lower 
than  computed.) 

(5)  First-order  cpproximations  are  used  to  describe 
the  operation  of  the  interferometer. 

For  both  the  single-  and  multiple-grating  ihterfero- 
metric  systems,  except  for  the  change  in  coordinate 
position,  we  have 

/(*)  ■  J*  ( f(a)da  +  J*  ^  /(a)  co*(biaMa, 

where  <*>  *  4 xFz/f. 

Consider  specifically  the  multiple-grating  processor 
as  illustrated  in  Fig.  10.  The  fringe  box  translates  down 
or  up  depending  on  the  position  of  the  input  point  For 
a  point  input  at  a  m  ao,  the  fringe  box  extends  along  the 
z  axis  up  to 

where 


T  •  un  ^in"1  |\F  -  tin  |un_l  jyj  (23) 

corresponding  to  a  cutoff  spatial  frequency  of 

,  2  FIA  2Dcto\ 

Mao)--(--— )•  (24) 

For  the  on-axis  point  «o  *  0,  the  cutoff  frequency  is 
a  maximum  and  is  given  by 


where 


0  ■  (261 

is  the  diffraction  angle  shown  in  Fig.  10.  Notice  that 
the  cutoff  frequency  is  wavelength  dependent  We  also 
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Fi*.  Jl.  Space-bandwidth  products  of  a  single-grating  interfero¬ 
metric  processor,  a  multiple-grating  interferometric  processor,  and 
a  coherent  optical  p*ocessor. 


input  length.  To  maximize  further  the  SBWP  of  the 
multiple-grating  interferometer,  it  is  necessary  to 
minimize  the  length  of  the  interferometer  2D.  (This 
is  equivalent  to  moving  the  input  closer  to  the  Fourier 
transform  lens  ic  the  case  of  a  coherent  optical  proces¬ 
sor.)  As  we  shorten  the  interferometer,  the  gratings 
eventually  collapse  onto  e  single  plane,  and  we  essen¬ 
tially  have  the  single-grating  interferometer.  (This  is 
similar  to  putting  the  input  right  in  front  of  the  trans¬ 
form  lens  of  the  coherent  optical  processor.) 

The  above  space-bandwidth  product  is  computed 
from  geometric  considerations.  A  second  factor  that 
can  limit  the  space-bandwidth  product  of  the  system 
is  the  dynamic  range  of  the  output  detector  array.  The 
bias  level  is  determined  by  the  integrated  input  signal 
level 

fl  ■  J*4/(a)da. 


point  out  that,  except  for  the  cosine  factor,  Eq.  (25) 
describes  the  highest  frequency  that  can  be  coherently 
processed  by  a  Fourier  transform  lens  with  aperture 
width  A,  focal  length  f,  and  the  input  placed  a  focal 
length  in  front  of  the  lens.  In  fact,  an  expression 
equivalent  to  that  of  Eq.  (24)  for  such  a  coherent  pro¬ 
cessor  is  given  in  the  Appendix. 

Now  consider  an  input  signal  of  length  21,  extending 
from  a  »  —l  to  a  ■*  l.  The  end  point  /  establishes  a 
cutoff  frequency  vc(l).  Although  one  could  define  a 
space-bandwidth  product  (SBWP)  given  by  the  prod¬ 
uct  of  the  line  length  by  the  maximum  cutoff,  i.e., 
4/MO),  it  is  more  meaningful  for  arbitrary  inputs  to 
establish  the  SBWP  only  as  specified  over  the  uniform 
(nontapering)  part  of  the  spectrum  as  given  by  the 
lowest  cutoff  frequency  vc(l).  Thus  we  define  the 
SBWP  as 


SBWP(l)  -  4l*f  (I). 

which  with  substitution  from  Eq.  (24)  becomes 
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(27) 


(28) 


uations  (24)— (28)  also  hold  for  the  single-grating 
rferometer  by  letting  D  »  0.  The  single-grating 
fferometer  can  be  considered  to  be  the  limit  of  the 
iple-grating  interferometer  as  D  approaches  zero. 
|wen  from  Eq.  (28)  that  the  SBWP  is  maximized  for 
0,  i.e..  for  the  single-grating  interferometer, 
is  worthwhile  to  compare  the  above  results  with 
for  a  classical  Fourier  transform  lens  used  as  a 
tent  optical  processor  (see  Appendix).  The  SBWP, 
r  *d  in  a  similar  way  as  for  Eq.  (28),  for  the  coherent 
:sssor  is 


coh(SBWP(2i)|  -  2//A-<*in|Un-'(M  -  2/I//II).  (29) 

wherje  the  input  extends  from  -A/2  to  21  -  A/2. 

In|Fig.  ll,  the  SBWP  is  plotted  as  a  function  of  the 
length  of  the  input  for  the  coherent  optical  processor 
and  for  the  two-grating  interferometer  geometries. 
Notei  that  for  the  multiple-grating  interferometer  and 
the  coherent  optical  processor,  there  is  an  optimum 


As/  increases  so  does  the  bias  leveL  Thus  the  number 
of  Input  points  that  can  be  processed  for  the  worst  case 
is  equal  to  the  dynamic  range  of  the  output  detector. 
The  space-bandwidth  product,  however,  can  be  much 
larger  if  the  transform  of  the  input  consists  of  only  a 
small  number  of  points  (e.g.,  a  sinusoidal  input).  In 
that  case,  the  space-bandwidth  product  of  the  system 
is  limited  by  the  system  geometry. 

We  should  note  that  the  same  dynamic  range  limi¬ 
tation  also  applies  to  a  coherent  optical  processor  if  the 
complex  output  is  to  be  read  out  with  a  detector  array. 
However,  there  is  more  flexibility  with  the  coherent 
system.  For  example,  if  the  input  data  are  on  a  carrier 
and  the  information  around  dc  is  unimportant,  the 
reference  beam  can  be  lowered  to  a  level  corresponding 
to  the  highest  signal  level  in  the  desired  portion  of  the 
output  spectrum.  With  the  interferometric  processor, 
on  the  other  hand,  the  bias  is  fixed  at  a  level  corre¬ 
sponding  to  the  dc  output. 

We  emphasize  that  unlike  many  other  optical  systems 
that  utilize  gratings  (e.g.,  the  achromatic  Fourier 
transform  system),  spurious  diffraction  orders  do  not 
necessarily  pose  a  problem  for  the  grating  interfero¬ 
metric  processor.  With  the  use  of  a  beam  splitter  to 
decode  the  output,  the  spurious  orders  have  minimal 
effect,  since  they  propagate  in  directions  away  from  the 
output  plane  as  long  as  the  angular  subtend  of  the  input 
is  less  than  the  diffraction  angle  of  the  grating.  This  is 
true  for  both  the  single-grating  and  multiple-grating 
arrangements  utilizing  beam  splitters.  For  a  multi¬ 
ple-grating  system  that  employs  a  tilted  grating  to  de¬ 
modulate  the  output,  the  spurious  diffraction  orders  are 
also  not  a  serious  problem.  First,  the  higher  diffraction 
orders  propagate  away  from  the  output  plane,  and  for 
the  right  arrangement  they  will  not  pose  problem.  The 
zero  order  does  propagate  toward  the  output  plane,  but 
it  can  be  suppressed  significantly.  For  example,  if, a 
thin  sinusoidal  phase  grating  is  used,  the  powers  in  the 
various  diffraction  orders  follow  the  Bessel  functions. 
By  adjusting  the  exposure  in  the  fabrication  of  the 
grating  to  where  the  zeroth -order  Bessel  function  goes 
to  zero,  the  dc  term  can  be  substantially  suppressed. 
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Furthermore,  the  zeroth  order  is  incoherent  with  the 
first  orders  that  are  being  combined  (large  path  differ¬ 
ences).  Thus  the  effect  of  the  dc  term  that  remains  is 
to  add  a  small  uniform  bias  to  the  output.  Spurious 
orders,  however,  do  create  serious  problems  when  using 
a  tilted  grating  to  decode  the  output  of  a  single  grating 
system,  and  such  an  arrangement  is  generally  not  de¬ 
sirable. 

VM.  Applications 

The  interferometric  optical  processor  can  be  config¬ 
ured  to  perform  operations  currently  accomplished  by 
other  optical  processors.  In  addition,  it  can  perform 
operations  that  cannot  be  easily  done  otherwise,  since 
it  operates  on  incoherent  light. 

For  example,  to  process  a  temporal  signal,  the  scan¬ 
ning  electron  beam  of  a  CRT  can  be  used  to  convert  the 
temporal  electrical  signal  into  a  spatial  formatted  light 
signal.  The  processed  output  is  received  by  a  detector 
array  which  is  read  out  after  integrating  over  the  time 
of  a  complete  scan  of  the  electron  beam.  Such  a  system 
can  operate  over  a  large  range  of  frequencies  by  varying 
the  scan  rate  of  the  electron  beam.  The  ease  in  ob¬ 
taining  the  real  and  imaginary  parts  of  the  complex 
Fourier  transform  output  gives  the  interferometric 
processor  a  potential  advantage  over  optical  systems 
using  AO  modulators  and  SAW  devices. 

The  interferometer  is  a  parallel  processing  device,  and 
its  capability  is  not  fully  utilized  when  used  to  process 
single-channel  temporal  signals.  For  applications 
where  parallel  inputs  are  available  (e.g.,  from  an  an¬ 
tenna  array),  the  inputs  can  be  used  to  drive  individual 
LEDs  in  a  linear  array.  Together  with  an  array  of  in¬ 
dividually  addressable  detectors,  the  system  can  be  used 
to  perform  n- point  discrete  Fourier  transforms  at  a  rate 
determined  by  the  frequency  response  limit  of  the 
LEDs.  As  compared  to  the  vector-matrix  multiplier, 
the  interferometric  processor  has  the  advantage  of 
simplicity  requiring  only  simple  gratings  and  coilima- 
tion  optics  instead  of  complicated  2-D  masks  and 
beam-forming  optics. 

More  important,  however,  is  the  ability  of  the  inter¬ 
ferometric  optical  processor  to  operate  on  completely 
incoherent  inputs.  The  system,  therefore,  has  the  po¬ 
tential  to  perform  Fourier  transformations  on  signals 
directly  from  the  real  world  using  light  emanating  from 
the  object  scene  itself. 

VIIL  Summary 

We  have  introduced  an  optical  processor  which  is 
capable  of  performing  complex  spatial  1-D  Fourier 
transforms  on  incoherent  optical  inputs.  One  possible 
implementation  of  the  processing  approach  using  a 
grating  interferometer  was  presented,  and  the  operation 
of  an  interferometric  optical  processor  is  experimentally 
demonstrated.  The  potential  ability  of  the  system  to 
directly  process  images  using  the  light  emanated  from 
an  object  scene  removes  limitations  imposed  by  input 
devices  and  opens  new  avenues  for  optical  processing. 
The  research  on  this  type  of  interferometric  processor 
is  still  in  it  initial  stages,  and  much  has  yet  to  be  devei- 
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Fig.  12.  Geometry  of  a  coherent  optical  proceteor  for  the  processing 
of  real  inputs. 


oped.  The  optical  configurations  presented  are  only 
examples  of  many  other  possible  implementations  of  the 
interferometric  optical  processor.  Alternative  system 
configurations  with  different  performance  character¬ 
istics  that  are  optimized  for  specific  applications  can  be 
developed. 


The  work  was  supported  by  the  U.S.  Army  Research 
Office.  Parts  of  this  paper  were  presented  at  the  In¬ 
ternational  Optical  Computing  Conference,  Apr. 
1983.22 


Appendix:  Space-Bandwidth  Product  o<  a  Coherent 
Optical  Processor 

For  an  optical  Fourier  transform  lens  to  process  the 
input  signal  without  loss  of  information,  all  the  light 
emanating  from  the  input  aperture  must  be  passed  by 
the  lens  aperture,  i.e.,  no  vignetting.  Consider  the  1-D 
optical  processor  illustrated  in  Fig.  12.  At  position  a 
the  cutoff  frequency  defined  by  geometry  consider¬ 
ations  is  given  by 

>r (a)  "  A-1  sin|Un~'M  “  2|a|)/2/!l.  (All 

Thus  the  cutoff  frequency  varies  along  the  input  aper¬ 
ture.  It  is  possible  to  define  a  maximum  space-band¬ 
width  product  via 

SBWPmn  »  4/i>c(0),  (A2) 

where  the  input  is  assumed  to  extend  from  —l  to  +f  and 
the  extra  2  comes  from  the  two  sidebands.  However, 
to  insure  that  the  input  is  processed  without  loss  of  in¬ 
formation,  the  highest  spatial  frequency  of  the  input 
must  be  limited  to  *.  (/).  Then  the  SBWP  becomes 

SBWPU)  *  iiH'U).  i  Alt) 

If  the  input  is  real,  the  transform  is  Hermitian.  One 
sideband  of  the  output  spectrum  can  be  discarded 
without  loss  of  information.  By  placing  the  input 
asymmetrically  from  -A/2  to  2 1-  A/2,  the  (one-sided) 
cutoff  frequency  becomes 

■  *inltan*'[^  “  20/ffl.  iA4i 

where  0  <  /  <  A/2,  The  SBWP  becomes 

^BWP'l/J  •  2/j»r(2/t.  i  A"> i 
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To  determine  the  optimum  input  size  for  a  given 
value  of  A  and  f,  we  let  21  m  nA  and  take  the  partial 
derivative  of  Eq.  (A3)  to  get 

<«> 

u  -  it*  <p 

and  of  Eq.  (A5)  to  get 

<*71 

(1  -  r»)»  P 

For  lenses  with  //Noe.  of  //2  or  slower,  n  is  nearly  0.5  for 
both  cases.  That  is,  the  optimum  SBWP  tends  to  be 
obtained  when  the  input  aperture  size  Ls  one-half  of  the 
lens  aperture  size. 
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